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ABSTRACT 


A computer program was developed to analyze 
Hovmoller Diagram data. 

The program uses Fourier Analysis to extract 
wave information from the individual reports. Fourier 
Anakysiswistusedmagaineto get antormationdas stownow 
these waves change with time. A 30-day forecast of 
the data is then made for the period covered by the 
next 60BReDares. 

The program was used on two sets of data 
obtained from the Atmospheric Environment Service, 
Department of Forestry and the Environment. 

The results indicate a possible relation 
between the means of the reports and the solar radia- 
tion. Some forecast skill was shown when comparisons 
were made to climatology and persistence forecasts. 
However, the error in the best forecasts was rather 


large for operational applications at a Weather Office. 
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In 1949, E. Hovmoller (1949) introduced the 
Trough-and-Ridge Diagram. Such diagrams show how the 
500-mb heights change with longitude and time. Oscil- 
lations or waves in these heights produce ridges and 
troughs wnicn can be-tracked for Several days, weeks, 
and in the case of long waves, even months. Forecasters 
at the Atmospheric Environment Service (AES) have used 
this diagram as an aid to forecasting the movement of 
troughs. ridaqes. and: energy -pulses.» In current: prac- 
tice, the forecast method is partly subjective and only 
indicates the possible locations at which rapid deep- 
eningeoresumracescyClOnNeS may OCCU 

A second method using Fourier Analysis has 
also been developed at AES, a numerical forecast scheme 
which predicts the 500-mb heights and, in this way, 
indicates the location and magnitude of possible deep- 
SHindeese Ata themtimesote Writind. ola tt lem1s knowmeabout 
the accuracy of this method (see Balshaw(1975)). 

This report examines two Fourier Analysis models 
for forecasting the 500-mb heights displayed in a Hovmoller 
Diagram. 


The first model assumes that the mean of the 
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heights and the amplitude of the waves are constant. The 
second model assumes that the mean and wave amplitudes 
“anvewitietimes § FOrecaSusmor UDm=vOmUnIniy days Were 
produced using the second model and Statistics are ca |- 


culated to show the “skill” of the forecast method. 
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Phase-angle versus Time for Wave Number 2 
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The Yearly Trend of the Southern Data Set 
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Plots of Amplitude Versus Time for Two 
Space-waves. (Northern Data Set) 
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INTRODUCTION 


Analysis of Data 

Every twelve hours, at 00 GMT and 12 GMT, two 
sets of seventy-two Hovmoller Diagram heights are pro- 
duced at the Canadian Meteorological Centre (CMC). One 
set of heights refers to a latitude belt from 50°N to 
75°N, the Northern Data set. The other set of heights 
refers to a latitude belt from 30°N to 65°N, the Southern 


Data set. 


i Sener glt oC UR ErOduatl on 

Each of the seventy-two numbers of the height 
set is obtained by taking a six-point average, along a 
selected line of longitude, of the 500-mb heights within 
CiemLatnLude wel te eine ines Otml Ong vUdesse lected rare 
fives ongi tude degrees apart. (ines first number i neche 
height set refers to 30°w longitude. “Each succeeding 
number refers to the longitude line five longitude degrees 
fomthe castmOfmetie: OreViOUSM line: mln und Suwa Ye, ey dpistelevere 
represents a set of heights averaged over a complete 


Zonaleband around the Earth. Figure: 1.1] shows =the sposa- 
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tion of the numbers and an estimate of the altitude of 
the surface at those longitudes for both Northern and 
Southern Data sets. The heights, reported by CMC, are 
in decameters (dam) with the hundreds digit deleted. 
Pia smmeans enliate eit) one: of Sthamnumbers 1s 50... then the 
actual average height’ for that line of longitude is 556 
decameters. However, for the Northern Data set, this lack 
of a hundreds digit could produce possible ambiguities 
Since a height coded as 99 could mean a height of 499 
decameters or 599 decameters. These 499 dam or lower 
heights usually occur in the Northern Data set during 
winter. This possible problem was removed by subtracting 
500 decameters from the actual heights. This changed 
some values to negative numbers. The height-sets were 
then punched on computer data cards and separated into 
the two groups, Northern and Southern Data sets. 

There are two different ways to analyze the 
height data, either by graphical or by numerical tech- 


niques. 


ae Graphtcal sechniques 


Hovmoller(1949) was the first person to propose 
a method to analyze his Trough-Ridge diagram. His method 
ASedescriped wii Balshaw oterc 109 7D) To este lleebenn-gused 
BOddY | Sein ise COCHITQUCECOUNSTS Us. On Pp HOt i ndistnG ss 6 baron 


heights using longitude as abscissa and time as ordinate. 
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Then equi-height contours are drawn on the plot. ~ Inter- 
vals of 10 decameters were used by Hovmoller, but other 
researchers have used different contour intervals. For 
example, Bjerknes(1957) used an interval of 4 decameters 
and AES forecasters use 6 decameters. Once contours have 
been drawn, areas of minimum heights (LOWS) and areas of 
maximum heights (HIGHS) are shaded in. At a specific 
time, the position of a height minimum or height maximum 
is called a trough or ridge, respectively. Lines connect- 
ing the various time positions of a trough or ridge are 
drawn to show the movement of these features. 

Two types of waves can usually be diagnosed in 
these diagrams, long waves {wave no. 1 to 3) and major 
waves (wave no. 4 to 7). Long waves are usually slow 
moving or even stationary in the Hovmoller Diagrams. 

Major or secondary waves usually have trough or ridge 

Hines ss lODiIngsdt. Varving: andles co. tne apscaissa. To 
forecast the height pattern, these trough and ridge lines 
are extrapolated past the most recent data. Constructive 
wave interference is expected at the location and time 
GUEWhICh me bWOULPO Ud NullaLMeSmOmwebWO sYicoe elinessintersect. 
ifesUWwounddgem|iiess intersect. nein tiGiels eexpectedmdaid:, 
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tion at the surface to these lows and highs can be expected 


about tive longitude dedrees Co) tie east OF tnese Teatures. 


- 7» 3 eae 


7 7 : 7 a 
-_ aL a i; ~ 
»  dpveenl” eee ere iarelate bina . 
- | lois | 
VOORE, Dell sel fama: vit "Ye vie tie igh as r 
it avy 5Aint! ques MWe Al Ane, Peay odd ae - 
sVOTAAN TED PL TK tes SEH | be ie Cape |) Wpee: , oan ian 
punt | Al Tit a i hue ; ' ned ee! hide 
: ese J 
is » lu? ; ; | Ai swe a : wane ; 
J 7 fy 
sy, wa i J rey oA inte et 
) sea 
' j ; 
77 SP at ) ' 34 Age ' Bt) an 
aii e ‘ oy : Lig 
x i 
' bi Vr ! 7 niu he } 
; 4 
f ” d 0 if 
: 1 \ ior py 
7 4 
4 } rai u 
e , ; A t ( i] 
' a ;4 @ Ki , i i ol 
= @:t 7 7 
_ t ‘ dap - : fT, if? & te 
va ' a > : a — : 
7 eo J “ee | ig 7 4) nih ated aid 
aa 7 ; _ a_i 7 a 7 
a - vv : 
ee eee 
| —~“ i 


’ ~¢ a _ ; — 
aha aduaaka ‘i 


| = ee 
Pree od i “ey m ay, yd res 
a Zz” - “0 1% 348 a0 ia y : 


La ai ” 


7 


Tag 


1.3 Numerical. Techniques 

As mentioned previously, the height data at any 
given time are repetitive in space for every complete 
circuit around the Earth. Most numerical techniques 
take advantage of this fact and analyze the data using 
a Fourier Analysis. This produces values for the ampli- 
tude and phase of the dominant waves in the data. A 
dynamic equation such as Rossby's wave equation can then 
be used to forecast wave motion. The Rossby wave equa- 


f1ON 1s used aby Balshaw eb..al (hos 5). 


1.4 Initial Model 

hie initial attempp at forecasting the heights 
for the Hovmdller Diagram also used a Fourier Analysis 
to extract space-wave information from the height sets. 


The formulas used were those of Jenkins and Watts (1968): 


k 
S(b)) = AG {2x Am) cos{o(h ln) yeteB(m)sinD( lm) )} 
m=1 


+ A(k+1)cos(D(L,k+1)) 
era 
where D(L,m) = “sy (L-LO)m 


* Page 255 Hess(1959) or page 144 Petterssen(1956). 
Another dynamic equation is the Cressman's equation 
(page 160 Petterssen(1956)). 
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LO = spatial origin of wave reference 
system = 937 


E> = Spatial location of the height, value 

in the CMCereporteset: 

L 1 tor thestfirst number in the set, and 
ieenow Bheilastenumber ian the set. 


L 


N = 72 = the total number of values in 
the set 


AO = the mean height of the set 
2B(m) = sine component of the amplitude 
of the space-wave having wave 


number m 


m = wave number = number of space-waves 
required to completely circle the Earth 


k = a 1 = maximum number of waves 


having a non-zero sine component 


A(k+1) = amplitude of wave with wave number 
are 


‘The amplitudes for these waves in space were calculated 


by: 
- 4f 
m ) 7 S(CBjicos (i) Giasm)s) 
L=1 
1 N 
B(m) = v2 S12 )Sol Tiel) laa alae 
L=] 


This form of the amplitude equations was con- 
verted to the amplitude - phase form by the equations: 
k 
SAU yp Se (NO 2S ok AMET Vetere COAL sil jarltGie 
m=1 


+R(k+1)cos(D(L,m) ) 


where 2R(m) = 2¥v A neyo Bit my = amplitude 
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Q(m) = arctangent(-B(m)/A(m)) = phase. 


This conversion requires that: 
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- 
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R(m)cos(Q(m) ) 
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-R(m)sin(Q(m)). 


The Rossby wave equation was not used to fore- 
cast movement of the waves. Instead, a simple rate of 


phase-change equation was developed. 


vin) = AB08) = 3 avctan CHORAL = MEE 


The derivation of this equation (see Appendix A) requires 
three assumptions. These assumptions are: 


1) The mean of the data set at time zero is 
assumed equal to the mean of any other time, 


2) the amplitude of the space-waves, 2R(m), at 
time zero 1S assumed wconscantestor ali other 
times, 


and 3) the phase speed V(m) is constant with time. 


Thesrest.of this chapter 1Sea) critical examina 


tion of these assumptions. 
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for Southern Data Set versus time. It can be seen in 
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Enisetrgupe, that the mean, AQ of the height-see 1s not 
constant with time. The error produced by a forecast 
made using Assumption I was as much as five decameters 
after fifteen days. This model would be required to fore- 
Cast fOr as, long as fifleenedayss to, one month. with an 
eccUGdCyeOTea > hUmddln 6 PheYehOremndiine INO) On ety eedeca— 
meters for just the mean of the height-set was considered 
too large. 

Examination of Figure 1.2 reveals two types of 
time fluctuations in the means.  sihev first type was. a 
long-term trendand the other fluctuation type was smal] 
amplitude "time-waves". 

The trend was assumed to have a yearly period, 
and to Ae a function of time of the form 

Rane 


AQ(t) = D +H cos(s35) RCS Ln 
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where= t= has units of —haht—dayse== ine coefficients D, HH, 
and C were assumed to be constant with time. These 
coefficients were then determined using a least-square 
error fit of the means of the height-sets. The deriva- 
tion of the least-square error formula is shown in 
Appendix B. 

As expected, the calculated trend does fit the 
diagnostic means well. However, forecasts of the trend 
beyond the limits of the diagnostic means had a tendency 


to drift away from the actual means as the forecast time 
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increased. The error produced was in the range of five 
decameters for a ten-day forecast. This error was re- 
duced as the number of diagnostic means increased. 
Figure 1.3 shows the differences Y(t) between 
Unemactual meansmused as diagnostic idata: s(t) and the 


calculated trend AO(t), so that 


¥(t)) = S Ct) =A 0 ie 


If these differences were small, then they could have 
been neglected. However, some Of these time deviations 
were as large as four decameters, and errors of this 
magnitude in the mean could not be ignored. 

All in all, Assumption I was shown to be poor 


fore tne cCesiRed forecasts penivods 


1.6 Assumption II: Wave Amplitudes 

The second assumption is that the amplitude 
of each space-wave is constant with time. Figure 1.4 
shows how the amplitude, 2R(m), of wave number m = 1, 
Chances withetime in the ssoutnern Dawa Sets 0Lher geuapns 
of amplitude versus time can be found in Appendix € and 
DE 

As expected from work with Hovmoller Diagrams, 
theafargestedmplitudes (10 to i5udam) occur for space 


waves that have wave numbers one-to four. However, the 
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time fluctuations of amplitude for these space-waves 

were larger and more frequent than expected. For example, 
the amplitude of wave number 1 changed from 5 dam to 

15 dam-over a period of four days. 

For Space-waves with wave numbers five to 
twenty-five, the amplitudes decrease as wave numbers 
increase. After wave number 11, the maximum amplitude 
observed was less than 2 dam. For this reason, space- 
waves Of wave numbers greater than eleven were ignored. 

The time fluctuations observed in these space- 
wave amplitudes could have been produced by the Fourier 
Analysis calculation and might not nave indicated a true 
time fluctuation in the waves in the atmosphere. These 
fluctuations were fairly smooth and continuous over two 
to four days of data, which indicates that a true dynamic 
process produced them. However, these time fluctuations 


were not consistent with Assumption II. 


las Assumption LLia Wave Phase 

mhe Wasteassumpaion.wthat) cne change of phase- 
andlemts constant Witnetime.. was also checked. ihe 
calculated phase-angle for each space-wave was plotted 
against time (see Appendix E). Figure 1.5 and 1.6 
show corrected phase versus time plots for wave 


numbers 1 and 2. It can be seen in Figure 1.5 that 
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Time for Wave No. 1 
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Figure ls6. Phase-angle Versus 
Time for Wave No. 2 
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for most of the diagnostic period, wave number 1] jis 
quaststationary, with a phase-angle of 5 Oe (Oniqin=at 
150°E longitude). For eight days, however, a rapid 
retrogression was detected with a speed of 150 km/hr or 
about 81kts. (See Appendix F for method of calculating 
wave speed). This unexpected time fluctuation of phase- 
angle would be handled poorly by the initial model. 
Other graphs for the phase-angle versus time also show 
rapid fluctuations of phase-angle with time, which would 
Lidn eat eewavVer speeds 01' omit On CUNKtS 000d DOUL) 25) tO 
Soni ne 

The question again presents itself. Does the 
calculated fluctuation represent a true wave motion in 


the atmosphere? The phase-angle calculation 


Q(m) = arctan(-B(m)/A(m) ) 


TSMSCHSLUIVeEWNens Ene aGacio =b im) /Atm) 1S close ta zero 
but insensitive when the magnitude of the ratio 
-B(m)/A(m) is greater than four. This problem would not 


explain the slow 360° shift shown in Figure 1.5 but it 


might explain time fluctuations in the other space-waves. 


Other problems associated with the amplitude- 
phase representation of the space-waves were: 


1) The computer arctangent calculation only 
allows for phase-angles between +180 degrees. 
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When the phase-angle has a magnitude near 180 
degrees, «then=a slichtetbimesT | tctuetiom aight 
cause a phase-angle sign change. This problem 
can be solved with proper computer programmimg 
but this was not attempted in this study. 


and 2) Long-time averages over a 10 to 20 day period 
are required to produce smooth phase changes. 
This also was not attempted in this study. 


These graphs show that not only is the phase 
navds LOeWOrkKe With Dit @alsO sthemassumperoneo, constant 
pnase VElOCILY dS NOCecCOrReat. 

These studies of the three assumptions required 
to use the initial model showed the assumptions to be 
poor when actual data are used. Also, the various 
problems associated with the phase-angle calculation and 
interpretation suggest that adjustments to the model are 


needed for useful forecasts. 


1.8 Final Numerical Model 

The initial model had to be changed to a new 
form that would reproduce the Hovmoller Diagram heights 
with greater accuracy. The quantities AO, A(m) and B(m), 
that were assumed constant with time in the old model, 
were assumed to be periodic functions of time in the new 
model. The form of this time dependence was, for A(m) 


and B(m) 
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A(m) > A(t,m) 


ily, 


NWAVE 
RA(m) +25 UPACIEMCOS\E Ct. | OA (lem): } 
IT=1 


NWAVE 


B(m) + B(t,m) = RB(m) + 22 {FB(I,m)cos(E(t,1)+QB(I,m))} 


where 


And, 


AO > AO(t) 


where 


ed) 


NWAVE = maximum number of waves used to fore- 
cast amplitudes 


eS ACES We Ui 


RA(m) = mean of cosine amplitude = 2 for 
Space-wave with wave number m 


RB(m) = mean of sine amplitude + 2 for space- 
wave with wave number m 


2FA(I,m) = amplitude component of time-wave 
having wave number I. (cosine) 


2FB(I,m) = amplitude component of time-wave 
having wave number I. (sine) 


QA(I,m) phase angle of forecast time-wave 


having wave number I. (cosine) 


QB(I,m) = phase angle of forecast time-wave 
having wave number I. [sine) 
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G,H = amplitudes of trend component 

2K(I,m) = amplitude of forecast time-wave 
of differences Y(t) mentioned in 
Sec tn Chee. 5 


QO(I,m) = phase of forecast time-wave of 
daftferencesaiut). 


Note that the space-wave amplitudes do not have a trend 


removal or a trend forecast. 
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CHAPTER =2 
BROGRAN ShOMFOREOAST SHE7 Gi is 
Figure 2.1 shows a block diagram of the two 
computer programs used to forecast Hovmoller Diagram 


Heights. 


Coe PCOCwola lc eeata COMection 


The first program collects data and updates 
diagnostic fields. As a new set of heights for time 
lee US) ead, incor tenes computer. «da Space Fourier Analysis 
is performed on it. The mean height, AO(t*), of each 
report becomes the last element in a vector, Y(t), that 
Containss prevaous means Of diagnosticwaata Tor time t=i 
COmte-aty — | whew em COM DON eINUS @Otmiu eaeS Pd Cel ai p | hau d Gau 
tes, A(t.m) and Bit mee pecome the eleven) elemenusedn 
the last row of two storage matrices, AAMPL (t,m) and 
BAMPL (t,m), respectively, The first element in the rows 
of these matrices represents the component for wave number 
m=1, and the last element in the rows represents the 
component for wave number m=11. 

Next, a least-square time analysis is performed 
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Block Diagram of the Two 
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trend equation in the form previously given in section 
iO einem trendeiSsCalctldtedm Or muemclagnosti cada ta 
¥(t*) and is extended outside the diagnostic period for 
a thirty-day forecast period. The difference, between 
the real means in vector Y(t*) and the estimated trend 
its) eS eCOMDU Led sa Midis BOM ecm ie eCCiLO me 7 

Pies i2onrpereporVcewascmm Ussiudn 1 One ta Memty a cilen 
tnree blank- cards were read into the computer in its 
place. The program will then set the space amplitude 
components A(t',m) and B(t’.,m) equal to 99.99 for the 
eleven wave numbers, and the mean Y(t') is set equal 
COmz ey Ope hOretne  CalcuatTOnOieehemaey ia ti Ont MOmpl ie 
trend, Z(t'), the missing data are replaced by a local 
time average. This time average is taken over the six 
(Orme fewer) On-Zzero. Values of Y(t). three on each “side. 
Genicot) Oe Demc Ome LC mao 1 Cmed Osco Us pUderO fen e 
fincstepart of thes neighterorecasteprodram 1S shown an 
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PessecondepamtrOterneoe pro dranmrwem cdakesud 
time Fourier Analysis of the amplitudes in each column 
of the two matrices AAMPL and BAMPL, and it takes a time 
Fourier Analysis of the deviations in Z(t*). On each 


occasion that a bime Fourier Analysis is carried out, the 
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resultant time-waves are sorted according to the magnitude 
of the time-wave amplitudes. Large-amplitude time-waves 
are listed first and small-amplitude time-waves are listed 
Tas c. 

At this point, the program can perform in var- 
ious optional ways. To determine which of these options 
is used, five parameters must be specified (see Appendix 
I for a complete description of these parameters). The 
three most important of these parameters are NWAVE, PRINT 
and JWAVE. The first NWAVE of the sorted resultant waves 
are used to forecast the ampiitudes in the columns of 
AAMPL and BAMPL, and to forecast the deviations Z(t). 

The second parameter, PRINT, is the most impor- 
tant. Jt alTows= forecasts to be made using persistence. 
climatology, unsmoothed time-waves, and smoothed time- 
waves. A persistence forecast (PRINT = 4.0) just repeats 
EDCeMOS Crecente lagnostic rrepoOntscVers Cher Tullerorecast 
period* A climatology forecast (PRINT = +5.5) is equi- 
valent to NWAVE = 0. This means that only the average 
OPeeachmool Umi lis euSed HCOsmaKe vas LIMGssfOrLecdSt. © O0heatne 
amplitude in each column of AAMPL and BAMPL. It also 
Meansmenat themaVverades Otc (tj) sismuSsedeto storecast. athe 
deviations from the trend. For unsmoothed time-waves to 


be used, PRINT must have a value greater than 6.0. If 


‘J erOreaspersictence, forecast, only: the VeriiicCationerro= 
gram needs to be run (see Appendix J). 
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PRINT 18 less than 3.5, time-smoothed waves are used to 
forecast the columns of AAMPL and BAMPL, and the devia- 
tions in Z(t). This smoothing is performed on the time- 
waves before they are sorted according to magnitude. 

The third important parameter, JWAVE, determines 
the number of space-waves (number of columns of AAMPL and 
BAMPL) to be used to forecast the heights. This parameter 
aso cd Wenavera VauUe fromelueLO il. 

Once the number of space-waves, JWAVE, has been 
specified, the strongest NWAVE time-waves forecast the 
first JWAVE columns of AAMPL and BAMPL and the deviations 
frome: urend “LUt). @ [new hercntse tor @ Torecast time 
(t+1) are then reconstituted from the forecast trend, the 
deviation from the trend Z(t+1) and the forecast space- 


wave amplitudes in AAMPL (t+l,m) and BAMPL (t+1,m). 
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SMOOTHING 


The time variation of the frequency of the 
time-waves is one of the major problems with this fore- 
cast model. Therefore, instead of the wave being repre- 
sented by a single frequency, as in the model, it should 
be represented over a frequency band, with the band width 
depending on the wave. This causes an error in the fore- 
cast of the space-wave amplitudes in the columns of the 
matrices AAMPL and BAMPL. The forecast program assumes 
that the time-waves, diagnosed in the columns of AAMPL 
and BAMPL, have a fixed wavelength during the whole fore- 
Ca Site per od. 

To try and minimize the problem of wavelength 
variations during the forecast period, a time smoothing 
was performed in Program II. Three different smoothing 


options were tested. 


3.1 Smoother Number 2 
Lachapelle(1977) suggested one method of smoothing 
the waves. In this method the mean-square contribution 


for wave 1. AAU! )se0is calculated! bya tnems oOrmue ar 


24 


gig Gn tahee) phe Ta ! aus | 
Huh “s BAIVD Ei he | athens raat EG we’! > 
wie patted ty ty win Ve Weg (eawravtet ; 
ete Wo teen Hap pi rumegiog’t nlindty ‘teaser 
Cap he hn AEN ey awh Aare eas b yee, eT | 
Te ey ba eye Ad? an gna 
; ‘7 ny | ike af pct eyeis od te or 
tn. do sdecne /eutval it eft tae A, 
‘than he’ melt ‘ in ig ba fale hd pote al} en 
$F alae val ght Wat OTHER IVS RKTT 5 avon. Athy 
| »Aarvad) 205 

Mic icekgltn Pie Ma “abatnre aha ped we 
Qura' pile GHP >. beluga Rap ai ns prengh can) Hate 
Si Vogts tnted plc ewih ° cll apfee yh WP hanspiagg oa 
| ee eo 


Mm 
nt 


ramplitude/2™= amoimtude/el = 2. 


Thewaverage of AA(])*0Ver three =waves 15° taken to be the 
smoothed mean-square. The frequency "f" is weighted 
uSing the mean-square as weights. Also the phase "Q" 
is smoothed by taking the average of the sine "BR" and 
cosine "AR" components. The arctangent of the ratio 
BR/AR is then taken to represent the smooth phase. 

This method was used as one smoothing option. 
The only change was due to a different definition for 
the phase. Lachapelle(1977) used the form Rcos(X-Q) in 
iis ime. Spectralwandlysises sin etniswprognan: the 4 Orn 
Rcos(X+Q) was used. This required the arctangent of 
the ratio -BR/AR to be used for the smoothed phase. 

To switch between amplitude-phase and sine- 
cosine forms, use the formula: 

BCI) =—-R(1)sin(Qul)) ARG) Saeki Ie RCI 

BX eee ly Cos.0 l)a)e 
For wave number I, smooth mean-square = {AA(I-1) 

+ AA(I) + AA(1I+1)}/3 
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smooth phase = arctangent (-BR/AR) 
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This formulation appeared to produce smooth 
phase-angles, but forecasts made using these smoothed 
angles showed phase shift problems when compared to the 
diagnostic data. Two other smoothing techniques are 
available in the program. Smoother No. 1 and Smoother 
Noweo careebothesilightavaniablonsmoreomeoeneralo. cc. that 
try to couple the calculation of smooth phase and smooth 


amplitude. 


3.2 Smoother Number ] 

In this option the frequency and phase were 
both smoothed as in Smoother No. 2. However, the mean- 
Square was calculated from a formula similar to that used 
to add three vectors in two dimensions. The resultant 
vector magnitude represents the smooth amplitude and 


its direction represents the phase. For wave number I 
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as before, but the smooth mean-square = 2(AR® + BR). 


This linked the calculation of the smooth mean- 


sauare and the smooth phase by using AR and BR in both. 
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smooth phase calculation and the smooth mean-square 
Calculation. | For the smooth phase calculation.) oa pro- 
blem arises since the mean-square is always positive and 
COMmeains NO phase ntormation. “slo counmmeract. thus 
problem the program uses the form 2(magnitude R)-R, in- 


Stead of the mean-square, to preserve the sign. 


For wave number 1: 
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smooth mean-square = (|AR] + |BR|)2 


smooth phase = arctan (-BR/AR) 


The following example demonstrates the results 


of these smoothing options: 


Initial Waves 


I Frequency f(1) AGL) B(1) RGTa) 0(1) AA(1) 
cycles cos Sin AMPL PHASE MEAN- SOQ 

a per hour _——__—s dam dam dam DEG DAM SQ 

] 0.00093 2 ns A eel Ae e026 Se elle 4.07 8 9206 
0.00185 -0.9484 042-38 260 49 WG Geom 200 el aoe 23 
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Result of the smoothing options, used for wave number 2 


ares 
F(1) A(T) B(1) R(1) Q(1) AA(1) 


unsmoothed 
0.00185 -0.9484 O22 30 2ee Uo oe LOD meg ee 


smoother #1 
OVOO0T34” "=0 35783 -0738983 07529 Pic eae 6 OU 


smoother #2 
0.00134 -0.61334 -0.66945 0.908 E3280) 0.64594 


smoother #3 
OT00T34 ~=0743896 =0763870 0475 Wea at el Oe eee 


sea Varlvances Equalization 

The total Variance is calculated from the sum 
of the mean-square for each frequency. This variance 
changes when the mean-squares are smoothed. The pro- 
gram attempts to make the variance after smoothing 
"SRMAX" equal to the variance before smoothing "SUM". 
This equalization is performed in the smoothing sub- 
routine by multipling the power fraction AA(1)/SRMAX 
by the variance "SUM" This adjustment technique does not 
bring the smooth amplitude to the same magnitude as that 


of the un-smoothed amplitude. 


3.5 General Comments on the Smoothing Options 


The greatest problem that occurred with the 
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smoothing method was the calculation of a smooth phase. 
As mentioned previously, Smoother No. 2 and No. 1 both 
appear to produce phase shifts of +20 degrees for the 
Strongest smooth waves relative to the strongest unsmo- 
othed waves. For a very strong wave, the amplitude would 
be reduced slightly by smoothing and the phase and fre- 
quency should not change. Smoother No. 3 produces little 
phase-change (+3 degrees) for these strong waves. 

Inmieach of, the smootningy optaonsmathe, first, fre- 
quency is smoothed by taking the average between the 
first and second frequency. The last frequency is never 
smoothed and is assumed so small that it 1S not even 
printed out (only the first 40 frequencies are printed 
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CHAPTER 4 


VERIFICATION 


Ae batho ele Se Used 

Sixty 12-hour forecasts were made for the 30- 
day period after the last day of diagnostic data. Several 
statistics were generated for each of these forecasts .* 


These Statistics wene:: 


We 
The mean of the forecast $0 = easy A) 
I=1 
1 die 
The mean of the real data YO =35- 2 Y(I) 
12 yn] 
1 TZ 
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Tie  Diascdmstandard sdevindllOneOreuie | 0recdst 


Thee biased Standard ceviadtion 0s) the real data 


sy= vq (y(1)-vo0)* 
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The mean absolute difference = v2 NY Gla) = Slo) 
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Here, the real data referred to were the actual 
CMC reported heights. These data were not used in the 
diagnostic data sets. Therefore, the statistics generated 
do not depend on the diagnostic data but show how the 
forecasts would verify in an operational forecast office. 
KDCMV al UCSHOMe ate eSita taS.b0/CS euuke wane O Pcie. 
the correlation coefficient and the mean-Square error of 
estimate, respectively, seemed to correspond well with 
subjectively determined good and bad forecasts. Since 
the forecast model was designed to produce "good" long 
range forecasts, the average R and the average MSEE of 
the first 28 forecasts were used as indications of 
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fe sles Sample 


To see if the whole program could forecast 
heights accurately, a test data set was produced. This 
test sample set contained 11 space-waves and each of 
these space-waves contained three time-waves.* A yearly 
trend was also added to the mean of each test report. 
The test sample repeated itself after 30 reports, except 
for the yearly trend. 

Near perfect forecasts were made when the pro- 
gram was used with diagnostic data sets of length 30 
reports and the time-waves unsmoothed. The length of 
the diagnostic sets, NUMB, was increased to 44 and 52 
reports and forecasts made. The results of these fore- 


casts are discussed below. 
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Figure 5.2 shows how the mean R and mean MSEE 
Stdatiusitiresmwene related Stor ethe iforecasts formvarious 
Sizes SU NOUMBY of the testediagnostic tdata ysets sikhe 


relation of these means was a linear one with 
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Since these two statistics were related, only one need 
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be discussed. For NUMB = 30, the maximum MSEE was 

Ae one over a nine-day period. This excellent result 
indicates that the program works for data which exactly 
Fit tne model; 

One of the problems of using a Fourier Analysis 
model is called wave truncation. This means that the 
wave being analyzed does not complete an integer number 
of cycles in the diagnostic data. This problem would 
cause small errors inethe forecastmat the truncated wave 
has a small amplitude. However, if the amplitude of the 
truncated wave is large and if the wave does not complete 
even one cycle then the error in the forecast can be 
large. 

To see how the program handles truncated waves, 
the test sample was used with NUMB = 44 and 52. For 
NUMB = 44 the longest waves completed only one and one 
half cycles. Table 4.1 shows how the values of MSEE 
changed with the forecast time for unsmoothed time-waves 
and smoothed (Smoother No. 2) time-waves. The values in 
this table show that smoothing the time-waves does reduce 
the error caused by wave truncation. 
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Porecast 


PROGRAM TEST USING SAMPLE DATA 


| NUMB = 44 NUMB=52 


Smoothed 
MSE 


Unsmoothed 


Time (% days) MSE 
1 42 
2 67 
3 106 
4 62 
5 32 
6 14 
7 38 
8 16 
9 53 

10 21 
11 73 
12 67 
13 92 
14 59 
15 82 
16 94 
i 126 
18 86 
19 62 
AVERAGE 
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CTT TerenimuUimess tne Values sin» lables4s teare dispaced 

So that they refer to the same times. For NUMB = 52, 
more data were available than for NUMB = 44. The MSEE 
values for these forecasts should be less than the values 
with NUMB = 44. This was not the case. Both sets of 
height forecasts used. smoother option No. 2, both were 
variance adjusted, and both had mean MSEE = 28 for the 
first 19 forecasts made. This indicates that the smooth- 
ing process needs more data, say several more complete 
cycles, before, a noticable increase ot forecast accuracy 
can be obtained. The change from 1.5 to 1.7 cycles 1s 


not enough. 


4.4 Other Smoothing Methods 


The smoothing options mentioned in Chapter 3 
were not the only smoothing techniques tested with the 


test sample data. Several other methods that were tested 


1) Smoother No. 2 with the variance not adjusted, 


2) Smoother No. 2 with the variance not adjusted 
and the mean-square of the first smooth wave 
found by 

AA(O) + AA(1) + AA(2) 
3 


where AA(0) = 0.00, always, 


3) Smoother No. 2 with variance adjusted but first 
wave not smoothed, 


4) Smoother No. 2 with all frequencies not weighted, 
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and 5) Smoother No. 2 with all frequencies not weighted 
and first time-wave not smoothed. 


The results of these smoothing methods are 
SHOWN neAppendinwt4s Teron thesemres toeum tetiss clear tnat 
the following produced lower average MSEE values: 
1) Variance equalization 
2) Weighted frequencies 
and 3) Smoothing first time-wave, using mean-square 


formula = aL Bees 
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CHAPTER 5 


RESULTS USING SS OUTHERNS OA Ties |) 
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Part 1 of the program was used with various 
values of NUMB (i.e. diagnostic data sets of various 
lengths) to test the trend calculation. The trend 


equation tested was; 


AO(t) = D+ H conte 


In lable 5-1, the values of. the coefficients (D, H and 
C) are shown for each value of NUMB. It can be seen 
that the value of these coefficients approach limits as 
the value of NUMB increases. To see how stable these 
limits were for different samples of the Southern Data 
set, tests were made with NUMB kept constant at 256, 

but the starting and ending time of the diagnostic data 
AdJUStECdsE Na DVel oO. Z2 ae SHOWSe the reSuUltSe nO eC) Gitano f 
these tests. The average value of the coefficients and 
the unbiased estimate of the standard deviation are 


shown below. 
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fap leea: Coefficients of Trend Equation. Mean 
Deviation from Trend, and Unbiased 


Variance for Various Values of NUMB 


D H C MEAN VARIANCE 
BAN dam dam dam dam eee 
5 1333 35/29 ne 0.000 02065 
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Table 5.2 Valuese fora lrends Coat ficients 
for Different Samples which 
have the same value of NUMB 


2) Ie SS Gs eo) fw ee Ce 
(Sal (02) tenk rl Cal (ey sey al 


Average Soo = Ome ie ale 


tes * 
ea saruaaiety 3 : m, my 


‘mye «a a ev. fee 403 aap 


40 


The maximum deviation from the mean for any of 
these coefficients was 2 dam. However, a steady change 
in value of the coefficients was observed as the diag- 
nostic data set was displaced in the direction of increas- 
ing time. Since the program only retains a maximum of 
300 reports for diagnostic purposes, some method of 
storage of trend information was required to carry out 
a more stable coefficient calculation. 

The method used in the program is to store the 
various parameters that are needed to calculate the 
COC Ti CLeEnTS, Dea Hoa and. Ce eWhenenew reports are read 
into the computer, the stored parameters are updated and 
then new coefficients are calculated. 

Figure 5.1 shows the trend values calculated, 
when NUMB = 298. The maximum, indicated by the trend, 
occured within two days of the actual maximum in the 
means and had a value within 1 dam of the actual maximum. 
The minimum value of the trend occured during the diag- 
nostic period and fit the means quite well. However, no 
statistics were calculated to indicate how well these 
trend values fit the means of all available reports. 

Also, a comparison can be made of the trend 
GULVYE.s( Figure 541), and) thencunve foOreunemnet solar, 
radiation arriving at the Earth's surface (see page 133 
in Hess(1959)). The two curves appear to be highly 


correlated, with an estimated time lag of 25 days in 
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The Yearly Trend, of the Southern Data Set 
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assumptions of a yearly trend. 


Two statistics were calculated to indicate the 


trendevalues sqoodnesss of fitetoutheameans of thesreports- 


These Statistics were: 


, NUMB 
the mean deviation from the trend MEAN = v2 Pits) 
j=1 
and the unbiased variance of the deviation around 
the trend NUMB F 
Variance “NUMB=T,", TEASE.) 


For this variance calculation it is assumed that the MEAN 
1s equal to zero since the least-square fit minimizes 
this mean deviation squared. Also note that all data are 
Sem O hue NCmeae CU lc bt On mm E mnOigSmal Net tear dudaese.t 
produced by estimating missing data report means, or 

any poor CMC analysis producing an abnormal mean are 


also used.* 
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Forecasts were made with the three smoothing 
options previously mentioned (NO. 1, NO. 2 and NO. 3) and 
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MSEE and average R were calculated. Since the use of 
two Statistics made it hard to compare the accuracy of 


fOrecCaS cS prasnew StatiStics wasaderined. 
me = Y okackaalor yi s NSEE- 


where r* has units of decameters-squared. 

Figure 5.2 shows the average R plotted against average 
MSEE OY Varlousetonecasts7=— he=ianess.go=through Rk = 
and MSEE = 0. Also, they go through the mean value for 
the forecast made, having the same NUMB. The statistic 
E~=nepeesents the distance, of the fTorecast, Values. of 
avegage Reand averages Moleg rome nm —-) 700 sand MSER =.) 700: 
Table 5.3 shows these values of r* for the _various fore- 
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As expected, the persistence forecasts give the 
DESte VESUlUSSOVeErEPer 1 00Sse0) | |ZacOmoOnnOUrnS. | HOWE Ver, 
persistence forecasts generally give poor results for 
DErHOdsS DEVOndesoehOlUrsee= ALSO. Trom Figure 5.2, 1t can 
be seen that the persistence forecasts have slightly 
higher correlation with the actual heights than would be 
predicted from the line through the means of forecasts. 
The persistence values also indicated how rapid the 
height field changed during the forecast period (large r* 


indicating rapid change). 
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Table 5.3. The Values of r* Determined From 
thewhn pst -2SeForecasts. ator 
Various Porecast0ntions 
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Where: NUMB = Size of diagnostic data set 
NWAVE = Number of time-waves used 
and JWAVE = Number of space-waves used. 
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Val viesmo tek andeMs BEatior chimatonogysforecasts 
indicate how rapidly the height field changes compared 
to the average over the diagnostic period. For a small 
diagnostic data set (NUMB = 90) the climatology results 
over the first twenty-eight forecasts had a lower r* 
thang ihesforecases usingethe largerediagnostic data set 
CNUMBE meee Oe alba Senseo Similar Chacac tend stice to: a 
persistence’ forecast Sincesthes "closer" diagnostic data 
set (NUMB = 90) produced the better forecasts. Since 
the climatology forecast for the group, NUMB = 240, was 
poor(i.e. r* large), then the forecasts with time added 
to climatology usually produced better forecasts (i.e. 
yee Sid 1 1). 

Looking at the R and MSEE values for the fore- 
casts using time-waves, it can be seen that unsmoothed 
time-waves gave good results for the initial forecasts 
but gave bad results for later forecasts, a consequence 
OF MnO tapering =tnemdiagnostic dataesetwat tne ends.) (On 
the other hand, the forecasts produced, using smoothed 
time-waves, had no preferred period for good or bad fore- 
casts. This was due to the change of phase, amplitude 
and frequency produced by the smoothing options. 


Usually the wave amplitude AAMPL (t,m) had one 
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or two well-separated strong waves in the data (for 

m = constant). The smoothing techniques, however, would 
cluster the smoothed waves (usually 3) around the fre- 
quency of each actual strong wave. When the program 
selected the three strongest time-waves, it usually 
got three smoothed time-waves that clustered around 
the strongest unsmoothed time-wave, and would not use 
any information from the second or third-strongest 
unsmoothed time-wave. The resultant wave amplitude 
produced from the combination of the three waves was 
larger than the original strong wave. The example 
below demonstrates this clustering in the smoothed 
waves. 

Example: These data are taken from the computer 
printout using actual height and wave- 
amplitude data. NUMB = 240 and the 
amplitude +2 being analyzed for time- 


Waves were ai) AAMPIBS((D as i ienencosiine 
component of space-wave number 3) 
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ae Mio (dam) (degrees) (dam) 
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Waves 500260 psitens 05 5 Sees) 
(Smoother 
Number 3) 


*k Frequency of wave with large amplitude, around which 
the smooth waves cluster 
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For this example, the maximum amplitude that these three 
Waves Can Obtain, if constructive interference occured, 
would be 1.6 dam for the actual wave and 2.1 for the 
smooth wave. This error in the amplitude in AAMPL, would 
produce a height error introduced not only by a phase- 
change, but also by an amplitude change in the space- 
wave. The program can't discriminate between the fre- 
QUencLes of @thersortedestrong@wave. « Lheretoresea tecan’t 
Skip from one set of clustered frequencies to another to 
USe -aS MUCK R IN foRmaki on ase ss avail ab) es 

Changes in the number of space-waves (JWAVE) 
used to produce the height forecasts appeared to produce 
only a slight change in accuracy of the forecasts (see 
bottom part of Table 5.3). For small JWAVE, an improvement 
WaSmesieen in the forecasts e9AISO. as. thesdiagnosticadata 
set size (NUMB) increased, a slight improvement was noted 
Ii sted WOreceSusee Hililse | Senndicacedan per CuTrem on eeby 
the grouping of the forecasts closer to a perfect fore- 
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CHAPTER 6 


RESULTS USING THE NORTHERN DATA SET 


G2 Sivrende Caden ations 

thes trend coeftacaents Dee and C4ewere poorly 
defined for the Northern Data Set. Table 6.1 shows that 
tne Values ofsthe coeffievencs continued to fluctuate 
even when the diagnostic data set was large. This lack 
Of FStability was due to thes large’ and rapid fluctuations 
in the means of the height sets (See Figure 6.1). The 
diagnostic data set covered the period Dec 1, O00Z, 1976, 
COmroeclins Onesie? melo / veel Seale ceed ey Det Caert lniSmnlis 
winter and, as expected, the heights are consistently low. 
It also meant that no noticeable spring warming had 
OECUrred@asS yet.@ sihhe resultant trendewas quiteut lat. (but 
did indicate a minimum mean around day 80 (Feb 18). This 
rsesliqhtly later than thevdate ofsthe minamum in the 
Southern Data set, and supports the assumption that the 


trend 1s produced by Solar heating. 
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Table 6.1. Values of Trend Coefficients and Mean 
and Variance Around the Trend, for 
Various: Values of | NUMB. 


NUMB D H C MEAN VARIANCE 
(dam) (dam) (dam) 

5 Sore -47.5 PReBENS 0.000 02107 
10 TPEOe6 -91.8 93.6 0.000 OR 92 
15 =1232:9 bia aya ral 10 2e56 0.060 Des 
20 -907.1 O77) Sleaze 0.065 O79) 
(ee A905 0 23005 Ole) On 24 12544 
30 -196.0 CaS L620 Omi 3 1.466 
35 G7, fog 8) -4.8 02095 2.054 
40 ioe 4 ©) -109.4 -20.8 O2094 22405 
45 Ose =1600.2 -27.9 0.106 AVES 
50 ii/eat ~1 53.47, -26.8 07095 Pest 
60 Ae -151.0 =) 0 0.078 LeU ge 
70 41.6 =2025 a9 Oe045 DE oU?e 
80 hea 0) oe taleee), G20 35 Gee ou 
90 324 -11.4 6.0 0.045 Sy 84/2 

100 ZOO -8.1 1S 0.058 68025 
110 -4.1 230 19.4 O2083 Fete ol 
120 -0.2 1929 feed, -0.013 EF MLNS 
130 1) lifel gre: re -0.007 The O Les: 
140 NG | ZA Z0R 2 0.016 6.885 
50 -8.6 Cal mee): eA! O20 179 7.040 
160 -5.0 2404 P99 0.006 6.885 
ZO -4.8 Pome 18.8 07005 6.566 
180 ea Teno Aaya 3) -0.003 Ue ewes) 
190 Gite 14.7 i779 -0.011 sha ea All 
200 sin ls LSieee LSe3 -0.010 S20 G7 
Cal. 0 igo 10m 9.0 -0.009 Saris! 
220 Loreal sy 18) SEA 0030 12.964 
230 229 eee Oral 0.024 14.412 
240 24.8 -0.2 -1.5 O03 145247 


Northern Data Set 
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did not verify as well as did those of the Southern Data 
SetUse Nable 6e2.qives,themvalucseot we -fogesevernal fore- 
cast options. These values are at least twice as large 
as those for the Southern Data set. Smoother No. 3 
Showed "increased skill" compared to the Climatology and 
Persistence forecasts while the unsmoothed time-waves 
showed “reduceumsk1) 1 > 

Part of -the fonecast enromewas duesto a sping 
Dad LLCrvesnii| bs as Can. De Seen ineit gure: ore. | Neathougn 
near JLOOW longitude starts to break down and to drift 
eastward. However, Some of the error was due to a poor 
trend calculation. The -error in@the mean started at 
S dam anda increased to e4udam attereco torecasts, and, to 
15 dam by the end of the 60-report (30 day) forecast 
period. 

The amplitudes of the waves (see Figure 6.3) 
were also larger in the Northern Data set than in the 
Southern set. The space-waves with wave numbers 1 and 2 
were still the most important and had amplitudes as 
large as 20 dam. The wave amplitudes decrease as wave 
numbers increase so that the amplitudes are generally 
below 5 dam once the wave number is greater than five. 
For any constant wave number, the time fluctuations of 
amplitude indicates long time-waves that might be re- 
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Statausties foraVarnious Forecasts Made 


Using the Northern Data Set with 


NUMB = 240 and JWAVE = 8 
Smoother NWAVE R MSEE r* 
Persistence -- -- Oeeia3 40-9 143.8 
Climatology -- -- 0.491 eden eye 
NONE 2 0.482 B37) leek 14626 
3 0.447 166.0 Oe 
1 0.483 sles 1227 
3 3 0.480 Ose 2620 
5 0.488 £6.20 IAS aes 
7 OF 42 aa 15924 
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FP sapeetsy (ey u6he Plots of Amplitude Versus Time 
for Two Space-waves. 


Northern Data Set 


WAVE NUMBER 1 


0.0 25.0 50.0 75 «0 100.0 125.0 
TD: SBI te) 


WAVE NUMBER 2 


NH a op 
Ca Hn 
minstomir ld 


ay 


0.0 25 .0 Heke. aia 100.0 125.0 
Teli ante yO) 


oo 


ie | 
x a) ge eek Daal ee BAI ie oe i fs 
vy fen dat i, al Hat 

Fy if .¥ Vite a in j + 


- 


ih + ' : ; 
as; rarer a qt cp ditere ee 


ge 


GHAP TER 7 


CONCLUSION 
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Figure 7.1 shows the actual Hovmoller Diagram 
andeFigure 7ec shows thee forecast diagwan aor the period 
Ma ricihimee. ted SLORAD OM pepe ee Ou emus I Gathe so OuUne Rn 
Data set. The forecast was made with Smoother No. 3, 
NUMB = 240, JWAVE = 11, and NWAVE = 2. These diagrams 
are Strikinglyeditrerenc. je forecast diagram is quite 
smooth, with the forecast ridges and troughs showing some 
resemblance to the actual ridges and troughs. However, 
this resemblance is not good enough to be used to make 
operational weather forecasts. 

This forecast diagram does show some skill 


Since it indicates three interesting features. These 


are: 
Le amstatuoneny «wndge sal ong 140°W, 
ciyeayretrogdressangqyridge sa long 7h aoe 

and Syeeraatrougn that splits into two troughs 


along 140 E. 
As stated above, these features are slightly correlated 
with the actual height changes. 
One of the things that is not observed in 
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the small number of time-waves used (NWAVE = 2); however, 
the use of more time waves produced smaller average cor- 


KeroaLIOn coethicients. 


7.2 Future Studies 

Oracheskir 1s cunrently stuaying, the, results of 
converting the 500-mb height data that are presented in 
a Hovmoller Diagram. He first converts the data to 
pstandavdized form. ‘hat stoned oreporte witheelements 
Xi, a mean AO and variance a the converted standardized 


report has elements Xi - AO , a mean of zero, and a 
2 


O 
variance of 1.00. The resulting data are then presented 
in the usual Hovmoller type diagrams. 

Using this procedure, he has found that ridges 
and troughs appear more distinct. Also, a diurnal cycle 
is easily observed in the data throughout the year. 

The height forecast program can be used to 
eheck if the Standardization vof theydata produces better 
forecasts than the actual data. Better results are 
expected. for forecasts of the standardized data since 
some of the "noise" would be removed. Also, since no 
trend would be observed, the error produced by the mean 
calculation would be removed. However, to reconvert the 


standardized data to 500-mb heights, the mean and the 


variance would also have to be forecast. The net effect 
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on forecast accuracy would have to be determined by 


experiments. 
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The forecasts for the Southern Data set were 
good for NUMB = 240 and should get better as NUMB in- 
creases (maximum NUMB = 300). On the other hand, the 
forecasts for the Northern Data set were poor. To 
produce better forecasts a much larger data set is 
needed. From Figure 6.1, a long wave of period around 
80 - 90 days, can be observed. This indicates that a 
Data Set close to one year in length (NUMB = 730) would 
be able to resolve this wave. The program would have to 
besadjusted Lo=operate- on such a large: mass of data, but 
the results should be much better than those presented 
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APPENDIX A 


Derivation of the equation for the phase speed V(m) 


Assume that the wave can be described by: 


AC tUcos DU. Cia mn) teh O eSani Odie) 


At) =. Ro CoS. @ (ity) 
B(t) = -Ry sin (o(t)) 
assume Cae Os Vito) 
where to: 6 Verereco. bimemzero i. er oe = 0) 


i] 


therefore A(t) Rcos(o +Vt) t 


R,cos(¢,)cos(Vt) - Rysin(o,)sin(Vt) 


AM Gees lm dio ly. 


B(t) = -R sin($,)cos(Vt) - R,cos(o,)sin(Vt) 


therefore 


Gis) a= A,cos(Vt) + Bysin(Vt) 


B(t) = B cos(Vt) - A sin(Vt) 
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Rearrange 


therefore 


Similarly 


therefore 


A(t) A 


= ~0 cos(Vt) + sin(Vt) 
Bo Be 
Bilets) ee B, cos(Vt) - sin(Vt) 
Ao Ao 
GOS) ean a - Bee ip rg x no + Fo 
0 0 0 0 
GOs ts A(t JA, + B(t)Bo 
oeeo Ne 
sin(Vt) = B At) - A B(t) 
aay + B Bo 
V = é arctangent ( B A(t) - A,B(t) 
BB t) + A At) 
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Amplitude versus Time for space-waves of wave 
number 1 to 11, respectively. Time is indicated as days 
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(Southern Data Set) 
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WAVE SPEED CALCULATION 


Mean radius of Earth 6371km = ne ‘ 


from axis = 


let no. of degrees phase change per day 


and the wave no. 


circumference = 27R 


wavelength 


Rcos45. = 4520km 
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get distance travelled/day by 
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APPENDIX G 


LISTING AND OUTPUT OF PART 1 OF THE 


HEIGHT FORECAST PROGRAM 


=oLORES) Old etrend = paramerers 


-Temporary storage for 


new trend Parameters 


-Program parameters - see Appendix I 
also see Format 7 line 193 


-Temporary storage for means of reports 


sNew "k" reports’ to be 
-Temporary storage for 
-Temporary storage for 
-Trend values storage 
-Deviations from trend 
=O Ldeab AMP lama traiix 

-Old AAMPL matrix 

-Old means of reports 


-Temporary storage for 
Darameters 


needed as input 
needed for output 
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updated AAMPL matrix 


updated BAMPL matrix 
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Format of New Data 

Three cards™are= required: 

Cordis Sera tse lert Ome Gee neni Meomch tS Ueyee 
columns and date in next 8 columns 
(see lines 36 and 37) 

Carded = sHeldnusmeo sto Scanner arsit 72. columns 

and any other parameter in the next 

4 columns 

Cards 3: Heightse qo tos7cesinetinsts72 columns 


and U,V (zonal and south wind). 
USon ds Vera remo Cause dain at Micsmeon; OG ware 
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HEIGHT FORECAST FROGRAM 


DIMENSION BUF (400) 

LWIMENSION DATAC72) vAAMF(300%11) »BAMFC300%11) 
INTEGER DATE sDAYy YEAR 

COMMON FREQ(151)AC151)*BC151) sAACIS1) 
COMMON /NO2/Y (300) 

READCS» 1K 

K=NUMBER OF NEW REPORTS TO BE REAL IN 
FORMAT (I3) 

CHECK STORAGE DATA BLOCK SIZE 
REALICL Ly 3IKLAZ 

REAICI2»3)KLAL 

FORMAT CIS) 

IF (KLA2.NE«KLAL)GO TO 500 

READIN 13» 4) LEN» Js Jy IM» JY 

LEN=NUMBER OF MEANS OF REFORTS STORED 

JY¥e JIM» JDy JH= YEAR» MONTH TAY yHOUR OF LAST REFORT 
FORMAT (IS, 413) 

WRITE (698) SY» JM» JDly JH 

FORMAT(’ ’» TATE OF LAST LATA’ 413) 

IF (KLA2.NE.LEN)GO TO S00 

IF(LEN.EQ.0)GO TO 20 
READC13 963) (RUF CL) yL=1*LEN) 
READIC12932) (CAAMP (Jel) eL=dyd1) 9J=19LEND 
REAICI 1932) (CRAMP (Jel) sL=4e11)9J=19LEN) 
TTIM=LEN+K 

REAL IN NEW REFORTS 

IF K AMOUNT OF REPORTS NOT AVAILABLE GO TO 15 
ITIM =TOTAL NUMBER OF REFORTS (OLD +NEU) 
Udie2® T=isk 

READIN S 9 SvEND=15) (DATACL) vyL=1924) yJHRy DAY s IMONy YEAR 
FORMAT (24F 3.0412) 

REAM CSS) (ATACL) 9L=25¥48) sDATE 
FORMAT (24F 3.0914) 

READS +7) (CLATACL) 9h. =49972) 9UyV 

FORMAT (24F3.00F4.19F 4.1) 

FRINT OUT NEW DATA 

WRITE (6216) YEAR» JMON + DAY » JHR 

FORMAT(’ “y’DATE OF NEW DATA %*4¢1Xy¥I2)) 
0 18 N=193 

ML=(N~1)k24 

MA= 14M. 

MB=244ML 

WRITE (4917) CLATACL) yL=MAyMB) 

FORMAT(’ ’y24F4.0) 

L=72 

SUM=0.00 
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OK 


32 
62 
63 
OK 


XL=FLOAT CL.) 

NOZRO=0 

CHECK FOR MISSING DATA 
NOZRO=NUMBER OF HEIGHTS OF VALUE ZERO 
N10 SO J=1lek 

TF CDATACIJ)D.-EQ.0.00)NOZRO=NOZROF1 
SUM=SUMTTATAC J) 

CALCULATE MEANS 

ANOT=SUM/ XL 

IR=T+LEN 

BUF CIK)=ANOT 

NO 12 Jel 

REMOVE MEAN FROM REFORT HEIGHTS 
DATACJO=DATAC I) -ANOT 

QEL=5.0 

IFCNOZRO.GT.35960 TO 13 

TAKE FOURTER ANALYSIS OF DATA 
CALL SPECCOATAyL +s DEL) 

NO 31 J=ivil 

IF CNOZRO-.GT-1O ACI) =99.99 

IF C(NOZRO.GT,. 10) BR( J =99.99 

UF -DATE DATA BLOCKS 
AAMF CIN s J2eACI) 

BAMF CIR s JD=BCI) 

CONTINUE 

MA=START FOINTs ITIMN=END FOINT 
MA=1 

ITT=ITIM 

DeGhTIM Pe) eS 00U GOTO ay 

IF NEEDED» SHIFT START ANT END FOINTS 
MA=ITIMN-299 


ITT=300 

WRITE (6960) ITI 
RORMAT GS iS 
JJ=0 

DO 21 L=MAvITTIM 
JJeJIt1 


RENUMBER ANDI STORE MEANS OF REFORTS 

YC JJ) =BUF CL) 

CONTINUE 

WRITE (46961) (BUF (L) »L=MAvITTIM) 

FORMAT(’ ’911F 7.3) 

WRITE (494) ITIL» JHRy DAY» JMONs YEAR 
WRITE C7 *62)01TI 

WRITE (862) 1TI 

WRITE (4963) (BUF (CL) xL=MAyITIM) 

WRITE (7932) (COAMF (Jeb) yL=teid)»J=MAyITIM) 
WRITE (8932) (CBAMF (Jel) eL=1911)9J=MAvyITIM) 
FORMAT CALIF 7.2) 

FORMAT (IS) 

FORMAT (11F 7.3) 

GO ON TO NEXT STEF IN FROGRAM 
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CALL WAVE CLEN+K) 

STOF 

WRITE CS*700) 

WARNING *X K VALUE TOO LARGE 

FORMAT (IHis’ DATA END HAS BEEN REACHED) 
Cte) Ome. 

WRITE (S9S0O1) 

ERROR ** DATA BLOCKS OF DIFFERENT SIZES 
FORMAT(’ ‘9s “KLA2y»KLALyLEN IN CONFLICT’) 
Sh 

END 


CKKRKKK TARE SPECTRAL ANALYSIS OF SCN) cook 


SUBROUTINE SFEC(S»NyDEL) 
DIMENSION S$(300) »R(150) »Q(150) 
COMMON FREQ(151)sAC151)»BC151) rAACISL) 
DNATA TOFI/6.2831853/ 

K=N/2 

Lens 2 toe 

FORMATCIS) 

FORMAT(’ ‘“’y»11F7.3) 

V=FLOAT CL) \ 
W=FLOAT CN) 

NO 3 M=1yk 

SUM1=0.00 

SUMN2=0,.00 

U=FLOAT CM? 

CALCULATE AMFLITUNE/2 

0 4 I=L 9N 

X=FLOAT (CT? 

D=TORTKUKCX-V)/W 
SUM1=SUMI+S (C1) *COS(D) 
SUM2=SUM2+S CID KSINGID 
ACH)=SUMI/W 

BRCM) =SUM2/W 

SUM=0.00 

FIND MEAN OF SCN) ANDI STORE IN ACL) 
Id S Js1yN 

SUM=SUM+S (J) 

ACL) =SUM/U 

SUM=0,.00 

FINI) MEAN-SQUARE OF SCN) MEAN 
AACLI=FACLIKACL) 

SUM=SUM+AA CL) 

RCLI=ARSCACL)) 

nO 7 J=1lyk 

ZJ=FLOAT (J) 

FREQ(J)=ZJI/ (HEL XW) 
RCJ=SQRTCACJ¥KACI) + BCJOKBCJ)) 
IF (ACJ) -EQ+0-00) ACJ) =0,0000001 
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157 
158 
ae 
160 
161 
142 
163 
164 
145 
144 
167 
148 
169 
179 
171 
£72 
173 
174 
a7 
174 
177 
178 
a7? 
180 
181 
182 
183 
184 
135 
186 
187 
188 
189 
190 
171 
eo 2 
193 
194 
oS) 
194 
a? / 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 


C *k FIND PHASE 
QCIO=ATAN2(-BC I) 9AC I) K360./TOFI 
C xx FIN! MEAN-SQUARE 
AACJIERER CIR (IK, 
IFCJULEQ.K)AACIO=AAC I) 2, 
C XxX SUM=VARTANCE=TOTAL MEAN-SQUARE OF WAVES 
SUM=SUM+AA CJ) 
IF(J.GT.40)G60 TO 7 
7 CONTINUE 
FREQCL)=0.000 
IFCACL).EQ.0-00)A(L)I=0.0000001 
QCLDY=ATAN2S (-0.002A(L))K3460./TOF I 
RETURN 
ENTI 


XKXX CALCULATE TRENID FARAMETERS ANT COEFFICIENTS 
USING LEAST SQUARE FIT TO SINE~-COSINE CURVE 
(SEE SAPPENIIIX RD 
SUBROUTINE WAVE CLENsK) 
JI= STARTING FOINT OF DLAGONISTIC DATA 
LEN=LENGTH OF BLOCKS=NUM 
KLAST= NO OF FOINTS OF ORIGINAL DATA 
J=POINT WHERE TREND STARTS 
ITEND=NQ OF FOINTS+i THAT WERE USED TO CALCULATE 
THE STORED TREND FARAMETERS 
ILAST=NO OF FOINTS TO BE USED FOR THE NEW TREND 
CALCULATION 
QNIMENSION -¥ (300) 4Y0€400) #22Z(300) 
COMMON /NO2Z/Y 
COMMON /SNOS/YO% ZZ 
DATA FI/3.141593/ 
READCSs7) IWAVE »yROsL ey DIT VyNWAVE »NUMB ye RMS » JT 
**K RMS=SQUARE ROOT OF VARTANCE OF DTFFERENCES 
FROM TRENT 
*¥K JI=START OF TATA IN DATA BLOCKS =1 ALWAYS 
xx L= NUMBER OF FORECASTS TO RE MATIE 
7 FORMATCISsF8.39ISesFS.39TSxTSeF 7.32715) 
Z=(2.XF1I)/730- 
IF (CLEN.NE.»NUMBIWRITE (6999) 
99 FORMATC’ ’s “LEN NOT = NUMB?) 
C *x READ IN STORED TRENT FARAMETERS 
READCL 44). TeWeSeFyQ 
REATIC1s 45 .EvFoGeANy TEND 
44 FORMATCSE1S.@) 
45 FORMATC4E1S.89IS) 
ILAST=IEND+K-1 
KLAST=LENtL+K 
J=LEN 
LEN=LEN+K 
C ** LEN NOW REPRESENTS LENGTH OF ALL DATA 
NUMB=LEN 
WRITE (46944) Te Ws Soke Q 
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WRITE (Sy 47) 
FORMAT (CS 4 AX)! 2X9! X91 4K VN 9 1SXs “COSYOUISK > 
1‘SIN’) 
nO 3 I=LENIyILAST 
J=J+1 
X=FLOAT CT) *Z 
U=COS(X 
=SINCX) 
Aaa 
IF (A,.EQ.0.00)G0 To 3 
T=T+U 
WeW+t¥ 
S=eSt (UK) 
Pek CUKUD 
Q=Q+ (VY) 
E=E+(AxU) 
FP + (UKA) 
G=G+tA 
AN=AN41. 
WRITE (5945) SeXrAvUyV 
FORMAT(’ ‘y9I5*4E15.8) 
CONTINUE 
WRITE (2944) Ty WeSo FQ 
IST=ILASTH1 
WRITE (2*45.EvFyGrANy IST 
R=TXS—FoKW 
l= ( CT XE-PG) & (SKS—FKQ) —RK (SKE-PUKF ) ) 
1/( (SkS—FKQ) KC TRT-PXAN) RRR 
C=( (SKE-PXE ) -RXLD /(SXS-FXQ) 
H=(E-TAD-S#C) /F 
PRINT OUT TREND) COEFFICIENTS 
WRITE (S24) DeHeC 
FORMAT(C’ ‘vE13.5%’ +/9E13.5+’COS +/sE13.5» SIN’) 
ILAST=ILASTHL 
CALL TRENDDCLENsIlyB»sCyILAST YL) 
CALI. DLFF CNUMB»LENs JI» RMS) 
WRITE (1497) LWAVE » ROyL y DLV» NWAVE » NUMER Y RMS» JT 
STOF 
ENE! 


TRENT FORECAST *x%xx 


SUBROUTINE TRENTIMCLEN» ty ByCrKLASTy JT» IENDy TLAST sl) 
DIMENSION YO¢(400) 9220300) 

COMMON /NOB/YO9Z2 
DATA PI/3.141593/ 
Z£2=(2.KP1)/730. 
FIND START-FOINT OF TREND TO FIT DATA BLOCKS 
MA=ILAST-LEN-L+1 

ITAST=ENT FOINT OF TREND OT MATCH DATA BLOCK 
TAST=ILAST<L 

ILAST=END-FOINT OF FORECAST FERIOL 
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WRITE C49 3) hMAyTASTyILAST 


FORMAT(C’ “»’ TRENT FROM’ sIéy’ TO’sIé2’ EXTENDE’» 
PG iéy” FOINTS = SPACER "EVERY 107) 


J=O 

GENERATE TREND FOR DIAGNOSTIC FERTOR ANT 
FOR FORECAST PERTON 

NO 1 T=MA» TLAST 

Jeb 

X=FPLOAT CI UKZ 

YOC J =D48KXCOS CX) +OKSINCX) 
CONTINUE 

STORE ANI PRINT OUT TRENT 
WRITE C994) I9MA 

FORMAT C213) 

Wile yee VeTOCL)s T=ty J) 
WRITE C692) CVYOCT) y T=12I210) 
BURMeN See he 72 5) 

RETURN 

END 


CALCULATE DIFFERENCES  XXxXxX 


SUBROUTINE DIFF CNUMB,yLEN» JT + RMS) 
DIMENSION Y¥(300)9X(400) 920300) 
COMMON /NO2/Y 
COMMON /NOB/Xso 2 
N10 2 J#JIyLeN 
IF CY¥CJI.EQ.0,.000060 TO 3 
Z(J)=7 C940) 

GO TO 2 


SMOOTHING MISSIMG DATA 
THIS IS A SUBROUTINE IN FART 2 OF THE 
HEIGHT FORECAST FROGRAM 
SUM=0,.000 

Me? 

MM=M 

A 7 POINT SMOOTHER FOR MISSING DATA 
no 4 T=1+MM 
Tecelrleiteoet reo GOs Wis 
TPt.GTs (LEN-3) 960 10 3 
N=J-4+1 

IF CY(N).EQ.0.000)G0 TO S 
SUM=SUMTY CN) 

GO TO 4 

M=M~1 

CONTINUE 

IF (M-LE.0)G0 TO 7 
SUM=SUM/FLOAT (M) 
Z(J)=SUM-X CJ) 

GO TO 2 

Z(J)=0.000 
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CONTINUE 

CALCULATE MEAN AND VARTANCE 
SUM1=0.000 

SUM2=0.000 

KR=LEN-JI+1 

N10 8 J=JIyLEN 
SUM2=SUM24(Z (J) KZ05)) 
SUMNL=SUM1I+Z (J) 

SMEAN=SUML/K 

VAR=SUM2/ (K-1) 

WRITE (412) IIT »LEN 

FORMAT(’ ‘’s “DIFFERENCE BETWEEN 12HR MEAN ANID TREND’ 


Liem ome Leleaee Oh ately 


JO=LEN-NUMB+ 1 

WRITE CLOs9)NUMB 

FORMAT CIS) 

FRINT OUT DIFFERENCES 
WRITECLOvs11)¢Z05) 9 JI=J09LEN) 
WRITE CGe911)9¢0Z0 5)» J=JDeLEN) 
EGitie «weir 74 5) 

WRITE C610) SMEANs VAR 
CALCULATE FMS 

RMS=SQRT CVAR) 

FORMATO’ “96X9 “MEAN=/ »F10.394X+ “VARTIANCE=% »F10.3) 
RE TURN 

ENS 
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APPENDIX H 


Sie NGRORSPART 2.0 bet HeESPROGRAM 


Files Needed 


Number 
IN 3 -Program parameters - see Appendix I 
also see Format 6 line 3/7 
IN a -Trend values storage 
(from Part 1 of the program) 
IN 10 -Deviations from trend - storage 
(from Part, yor the program) 
IN ne -Updated BAMPL matrix 
(from Part 1 of the program) 
IN V2 -Updated AAMPL matrix 
(from Part iMofethe program) 
OUT 14 -Storage for forecasts of the AAMPL matrix 
OUT 15 -Storage for forecasts of the BAMPL matrix 
OUT 16 =Storage: fore TOrecastSmoustne deviations 
from the trend 
OUT ES) eStorage: tor forecasts Of heignts 
Where 
IN needed as input 


ou 


OUT needed for output 
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FART 2 kx 


DIMENSION BUF (400) 

COMMON /NO4/Z (400) » AAMF' (400911) » RAMP (400511) 
COMMON FREQ(151)eRC151)°Q(151) 2AACIS1) 

CHECK SIZE OF STORED DATA BLOCKS 
READCL 1» 1 NUME 

REATIICLOs 1 NUMBL 

REANCL 2:1 NUMB2 

FORMATCIS) 

IF (NUMB1.NE.NUMB2)60 TO 900 

IF CNUMB1.NE+NUMB)GO TO 900 

GO TO 2 

WRITE (4% 100) NUMB»NUMBI »NUMB2 

FORMAT C2 % “HATA INABLOERS Alby TI PREREND SS (Zev 3 io) 
STOF 


READ IN DIFFERENCES (MEAN-TRENT) 
READCLOs 3) (20 J) 9J=1»NUMB) 

FORMATC’ ‘911F7.3) 

READ IN STORED AMPLITUDES/2 

READCL 214) (CAAMP CJsL) eL=1211)*J=1»NUMB) 
READCL1:4) (CRAMP (Jeb) »L=1e11)»J=1sNUME) 

FORMAT CLIF7.2) 

READ IN FROGRAM FARAMETERS 

IWAVE=NUMBER OF SORTED STRONG WAVES TO BE PRINTED 
FRINT=OFTION INDICATOR 

L=NUMBER OF FORECASTS TO BE MADE 

NIV=JWAVE IN REAL FORMAT 

NWAVE=MAXIMUM NUMBER OF TIME WAVES IN FORECAST 
NUM=NUMB=LENGTH OF DIAGNOSTIC DATA BLOCKS 

RMS» JI=FARAMETERS USED IN FART 1 

REANC3Ss6) IWAVEsFRINTs Le DIVeNWAVE »NUM es RMS» JT 
FORMAT CISsF8.3xISrF8.39TSrISvF 7.3215) 

WRITE (469172 NUME 

FORMAT(’ ’»’NUMBER OF DATA POINTS =’9I5) 


KLAST=NUMBt+L 

AMFPLO=0.15 

AMFLO=MINIMNUM AMPLITUDE ALLOWED FOR A FORECAST WAVE 
DEL=12.0 

DEL=NUMBER OF HOURS BETWEEN REFORTS 

START TO OPERATE ON DIFFERENCES 

WRITE (G97) 

FORMAT(’ ’*’SFEC FROM MEAN-TREND DIFFERENCES *% ) 

CALL SFPEC(Z»NUMBy DEL) 

CALL SFOUTCNUMBy RO» JHs SUMyFRINT yRMAXs JJ) 
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IF CARS(CFRINT).GT.3.5)G0 TO 5 

CALL SFCSM(JJ»SUMsF RINT) 

CALL. SORT CIWAVE » JHsSUMsFRINT ) 

CALL FCSTCNWAVE »ROrL y DEL» JLT»NUMB» RUF ys AMPLOsFRINT) 
NO 18 T=1»KLAST 

Z£¢T)=BUF CI) 


START TO OFERATE ON AMPLITUDES 

TA=WAVE NUMBER COUNTER 

FIRST AAMFL 

TA=1 

JWAVE=TTV/1 

DNATMS=99 .99 

WRITE CG7110TA 

FORMAT(’ “9 ’SFEC FROM COS AMPLITUDES’ sSX»’WAVE NO’:I5) 
NO 10 I=1»sNUME 

BUF CI) =AAMF CI» TA) 

CALL SM?7FTCBUF »yNUME»y DATMS > 

CALL SPEC CBUF »yNUMBsTEL) 

CALL SFOUT CNUMB: RO»: JHs SUM» PRINT ss RMAXs JJ) 

IF CABS(CF RINT) .«GT.3.5)60 TO 8 

CALL SFCSM(CJJsSUMsFRINT) 

CALL SORTCIWAVE » JHySUMsFRINT) 

CALL FCSTCNWAVE »yROsvL ye TEL»y JT »NUMBy BUF yAMFLOe PRINT} 
NO 19 T=1lyKLAST 

AAMF (Ir TA)=RUF CT) 


NEXT BAMFL 

WRITE CG¥14)1A 

FORMAT(’ ’+’SFPEC FROM SIN AMPLITUDES’ »8X*s “WAVE NO’ IS) 
00 13 I=1»NUMB 

BUF (I )=BAMP CT» TA) 

CALL SM7FT  CBUF »y NUMB» DATMS > 

CALL SPEC CRUF »NUMEBs TEL) 

CALL SFOUT CNUMBs RO» JHy SUM» PRINT» RMAX» J) 
IFCABS(PRINT).GT.3-5)G0 TO 9 

CALL SFCSM(JJ»sSUM*sPRINT) 

CALL SORTCIWAVE» JH» SUM» FRINT) 

CALL FCSTCNWAVEsROsL + DEL» JT »NUMB» BUF y AMPLOv PRINT) 
NO 20 T=isKLAST 

BAMF (I 9 ITA)=BRUF CT) 

CHECK TO SEE IF MAXIMUM WAVE NUMBER REACHED 
TA=IAt1 

IFCIA.LE.JWAVE)GO TO 14 

FORMAT (IS) 

STORE NEW DATA BLOCKS 

WRITE (14925) KLAST 

WRITE (C1525) KLAST 

WRITE C16s2S5) KLAST 

No 24 IT=1LyKLAST 


ey. ay SUnnee ebb e USRUDTT USAT OOD EES: Coane! 


(AWE Hai NYRI 
(5) 9 ee oe 
(id. vi ptlils cnden adnan. ome Arie 
~ 97Q) Gott tap Pes 
rTAte ue lanes. 20 - 
4 art ae ey BES I I OE 
[THINS HN WS 5 Grice, Las eh a4 Ce BD FT 


“CEL (om) SEAS TEs WERT PU ETE mnie gat "TLAMAYS OE ; 
_ 


CTA « aon ee . 
CAs. ere ME DT) ZA, a. 


Clb Birt hea 


lt he 


)PATOR ue es ay AME TROY Ag 1 
CHESS Wee RON RAE Saw TRON AIT md 


sits: 


“atyt2 


VRasHelel #e ne 
(Lytham EAL SI4nde FY kas 


sai ea eh Ha 
arb pd) Wt he 7 


hie Fedo) Ot 
pare? Wirr@uc ewe e2 _ 


(Ae GMit) Hare tu, a 
Vath. (TALMVOMAY AT « / 
9. AAAI AD 


Tia {Aiea of 00 
‘erie 


ee 2 9 OS Ff 24 OS $H OF $4 SF SF OO SH SH BS SH SH SS SH HO OH FH SH HH BH SH OH SH SH OH $4 SH HH HH FH SH FH CH HH FH $4 $H BH HH FH FH $4 $H TH OCH $4 OH FF 


107 
108 
109 
110 
114 
112 
113 
114 
115 
116 
7 
118 
119 
120 
124 
122 
133 
124 
125 
126 
127 
128 
129 
130 
1314 
132 
133 
134 
135 
136 
17 
138 
139 
140 
1414 
142 
143 
144 
145 
146 
147 
148 
149 
150 
1Si 
152 
153 
154 
155 
156 
157 
158 
159 


132 


WRITE (14923) (AAMP (Ts J) yJ=1» WAVE) 
WRITE (15923) (BAMP (Ly J)» J=1* WAVE) 
WRITE (16%23)(Z(5) 9 J=1*KLAST) 

CHECK TO SEE IF WANT AAMFL PRINTED 
IFCABSCFPRINT).LT.9.5)9GO TO 26 

NO 21 IT=NUMBsKLAST 

WRITE (4923) (AAMP CI» J)» J=1s WAVE) 
FORMAT(’ ’911F7.3) 

CONTINUE 

CONTINUE 

FORECAST REFORTS 

CALL FCAT (JWAVE »sKLAST »NUMB sy PRINT) 
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USED TO ESTIMATE MISSING ITiATA 


SUBROUTINE SM7PT CBRUF »NUMByTIATMS) 
NUMB= SIZE OF BUF (LESS THAN 400 ELEMENTS) 
IATMS=VALUE ASSIGNEN TO MISSING IATA 
DIMENSION BUF (400) 

nO i J=1»NUME 

IF (ARS (BUF (J)-IIATHS).GT.0.01)GO TO 1 
SUMN=0.00 

Me7 

DO 1=197 

N=J-441 

IFC(N.LT«1)+OR. (N.GT.NUMB) GO TO 3 

IF (ARS (BUF (J)-IATMS) «LT.0.01)G0 TO 3 
SUM=SUM+BUF (N) 

GO TO 2 

M=M— 4 

CONTINUE 

IF(M,LE.0)GO TO 4 

RUF (J) =SUM/FLOAT (M) 

SINCE THE MISSING [IATA IS REPLACED WHEN 
DETECTED THEN IT IS USED FOR THE NEXT 
REFORT IF IT IS ALSO MISSING 

GO TO 1 

BUF (J) =0.00 

CONTINUE 

RETURN 

ENII 


CKKKOKK PRINT OUT DATA FROM SPECTRAL ANALYSIS xxx 


SUBROUTINE SPOUT CNUMEy RO» JHySUMs PRINT» RMAXs JJ) 
COMMON FREQCI51) ¢RC1S1)°QC151) sAACIS1) 
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XX NUMB=SIZE OF DTASNOSTIC RATA SETS 

RO=MEAN OF REPORT 

JH=MAXTMUM NUMBER OF SPECTRAL FREQUENCIES 

RMAX=SUM OF THE SFECTRAL AMFLITUIES/2 

SUM=SUM OF THE MEAN SQUARES AACT) 

Q=PHASE 

R=AMPLITUNE/2 

JJ=NUMB/2 +1 

JH=JI-1 

SUM=0.00 

RMAX=0,000 

NO 18 I=1»JH 

RMAX=RMAX+R CT) 
18 SUM=SUM+AACT) 

IFC JH.GT.40) JH=40 

TPCPRINT EI. 0.900G0 TUS 

WRITE (6914) 
16 FORMAT(’ %’910Xs’FREQ’ 10% “AMFPL’ ¥10X+ “PHASE ’ 8X» 

1’MEAN SQR’»8Xs“FPOWER’ > 

WRITE CS 99 FREQ CII) sR OU 2 QC J) FAACII) 

00 8 I=1+/JH 

FOWER=AACT)/SUM 
8 WRITECGs9) FREQC LT) oe ROT)» Q¢01) sAACT) sFOWER 
9 FORMAT(* %sSF15.5) 
3 CONTINUE 

RO=R CS) 

IF ((€180.00-ABS(Q( JJ) )).LT- (0-10) )RO=—-RO 

SRM ic ota MAD 
2 FORMAT ¢’ y’VARTANCE =’ 9F7.395X» “MAXIMUM AMPLITUDE “ 

Lez a= oe 0nd) 

RETURN 

END 


al tal wey tee) (oa) teal toa 
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C¥KK*XXK TAKE SFECTRAL ANALYSIS OF SCN)  Xx*KKK 

is 
SUBROUTINE SFECCS» Ne TEL) 

DIMENSION $¢400) sAC150)»BC150) 
COMMON FREQC151)eR¢151)9QC151) sAACLTSI) 

*K ACI)D=COSINE COMPONENT OF WAVE NUMBER I 
E(I)=SINE COMPONENT OF WAVE NUMBER I 
N=NUMBER OF ELEMENTS IN VECTOR S (N LESS THAN 400) 

*K NOTE *K SQUARE OF MEAN OF S(N) STORED IN AACL) 

VATA TOFI/46.2831853/ 

K=N/2 

L=N/2 + 1 

FORMAT (IS) 

FORMAT(C’ “911F7.3) 

V=FLOAT CL) 

W=FLOAT (CN) 

NO 3 M=1yk 

SUM1=0,00 

SUM2=0,00 

U=FLOAT (CM) 
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[0 4 I=19N 

X=FLOAT CI) 

D=TOPTXUKCX-Y)/W 

SUM1=SUM1I+S(T) xCOS(D) 
SUM2=SUM2+5 (1) *SINCH) 

CALCULATE SINE AND COSINE COMPONENTS 
ACM) =SUM1/W | 

BRCM) =SUM2/U) 

SUM=0.00 

NO S J=1sN 

SUM=SUMTS (CJ) 

ACLI=SUM/W 

SUM=0.00 

AACLI=ACLIKACL) 

SUM=SUMNtAACL) 

RCLI=ARSCACL)) 

NO 7 J=1sk 

CALCULATE THE FREQUENCIES IN CYCLES FER HOUR 
ZI=FLOAT (I) 

FREQ(J2=ZI/7 (HEL KW) 
RCJ=SQRTCACIJORAC IO + BCJIKBCI)) 

IF CAC ID -EQ20.00) AC JI=0,0000001 
CALCULATE PHASE 

QC JI=ATANZ(-BCI) se ACI) *%560./TOPT 
CALCULATE MEAN SQUARE 

AAC IJORR CJ RR CIO K2. 
IFCJ,EQ.KIAAC JO FAACI) 72. 
SUM=SUMPAA CJ) 

IF(J-GT.40)GO0 TO 7 

CONTINUE 

FREQCL)=0,.00 

IF (ACL) -EQ-0-000) ACL) =0.00000001 
GALCULATE STHE FHASE OF THE MEAN OF SiN) 
IF Q=+»-180. THEN MEAN=~MEAN 

IF Q=0. THEN MEAN HAS NO SIGN CHANGE 
Q(L)=ATAN2(-0.002A(L) ¥360./TOFT 
RETURN 

END 


C 
CKKKKKK LIST WAVES ACCORDING TO AMPLITUDE MAGNITUDE 


feaj ispi tealisah tee! Ga) eB 


KX 


LARGEST AMPLITUDE WAVE LISTED FIRST 
SMALLEST AMPLITUDE WAVE LISTED LAST 
SUBROUTINE SORTCIWAVE»yL»SUM»yPRINT) 
IWAVE=NO OF FREQUENCIES 
L=NUMBER OF FORECASTS TO BE MADE 
SUM=TOTAL OF MEAN SQUARES OF WAVES 
PRINT=FPROGRAM FARAMETER 
ITOF=ROW INDICATOR 
COMMON FREQ(151),AMPL(C151) 2FHASECIS1) sAACIS1) 
ITOF=1 
STORE FIRST FREQUENCY 
BUF=FREQ(ITOF) 
BUA=AMFL (ITOF) 
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BUP=FHASE (I TOF) 

BUFR=AACTITOF) 

TTOF=ITOF+4 

SEARCH THROUGH ROWS 

DO 1 J=ITOP sk 

IF CBUA.GE.AMPL(J)2GO TO 1 

IF AMPLITUNE/2 LARGER THAN STORED WAVE THEN STORE 
THIS WAVE INSTEAD AND MOVE OLD STORED WAVE TO THIS 
ROW 

STORF=FREQ( J) 

STORA=AMFL (J) 

STORF=PHASE (J) 

STOFR=AA CI) 

FREQ(J)=BUF 

AMPL (JO =BUA 

FHASE (J) =BUF 

AACJ)=BUFR 

BUF=STORF 

RUA=STORA 

BUF =STORF 

BUF R=STOFR 

CONTINUE 

K=ITOF~1. 

SHIFT STORED WAVE TO FREVIOUS ROW 
FREQ(K)=BUF 

ANFL CK) =BUA 

FHASE (RK) =BUF 

AA CK) =BUFR 

IF(CK.EQ.IWAVE GO TO 6 
IFCK.GT-LIGO TO S 

GO TO 4 

WRITE C&s7)K 

FORMAT(’ “»sI3%’ THERE AFFEARS TO RE SOME ERROR’ » 


hee Le Ba ee) 


RETURN 

CONTINUE 

CHECK TO SEE IF WANT STRONGEST WAVES FRINTED OUT 
IF FRINT GREATER THAN 0.5 THEN THEY ARE FRINTED OUT 
IF (CARS(PRINT)«LT.O.S)GO TO 3 

WRITE (699) TWAVE 

FORMAT(’ ’9% THE ‘9IS*’ STRONGEST WAVES”) 

WRITE (6913) 

FORMATO’ %9’ FREQ’ s3X*’AMPL’ 4X» “FHASE’ ¥6X% “RMS 7) 
CONTINUE 

TOT=0.00 

AL UF MEAN SQUARES OF THE 6 STRONGEST WAVES 

No 8 Jails TWAVE 

IF(J.GT.6)G0 TO 12 

TOT=TOTHAA CJS) 

CONTINUE 

IF (CARS(FRINT).LT.0O.5)GO TO 8 

WRITE (692) FREQ(J) sAMFPL (J) s PHASE CJ) rAACJI) 

FORMAT(’ “9F7.SrF 8+ 372k Pe 3) 
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8B CONTINUE 
TOTAL=TOT/SUM 
WRITE (6210) TOTAL 
10 FORMAT(’ ’+’FOR FIRST 6 WAVES VARIANCE EXPLAINED’ 


LF 7.3) 
RETURN 
ENII 
C 
CK FORECAST TIME WAVES IN AAMFLy BAMPL» 
: ANI DIFFERENCES 4% k4X 


SUBROUTINE FCSTCNWAVE »ROyL » TEL » JT »sNUMEy VALUE » 

LAMFPLOsFRINT) 

DIMENSTON VALUE (400) sACIS1) 

COMMON FREQ(C151) sAMFL (C151) sPHASECIS1) xAACISI) 
DATA FI/3.141593/ 

Q=2.XFI/340.- 

X=2,. KP TXTEL 

KLAST=NUMETL 

K=90 

C %’k BY-FASS THIS LOOF IF CLIMATOLOGY FORECAST WANTET! 

IF CABS (CARS (FRINTI-S.S).LT.0.1)960 TO 8 

NO 7 IT=isNWAVE 

ACT) =2.XAMPIL CT) 


K=K+1 
C. ** CHECK TO SEE IF THE AMPLITUDE IS LARGE. ENOUGH 
C TO BE USE $ IF AMFL LESS THAN AMFLO THEN NOT USED 
IF CAMPL(I>.LT-AMFLO)GO TO 8 
7 CONTINUE 


8 WRITECGs10)K 
10 FORMAT(’ ’»’ THE FIRST ‘’sISyr’WAVES WILL BE USED’) 
JORG=1+ (NUME/2 ) 
C *k JORG=ORIGIN 
NO 3 JeJIyKLAST 
SUM=RO 
C *xk CHECK FOR CLIMAT FORECAST 
IF CABS CARS CFPRINT)-~S-5).+LT.0.1)60 TO 6 
JN=J-JORG 
T=FLOAT CIN) 
C x#%k AIL FORECAST WAVES TO MEAN OF AMPLITUDE 
10 4 T=19k 
SUM=SUM+A CT) XCOS(XXFREQ( I) XT4+QXFHASE (CT) ) 
CONTINUE 
VALUE (J) =SUM 
CONTINUE 
RETURN 
END! 
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C 

CKKKKKK FORECAST REPORTS  **xKKX 

C 
SUBROUTINE FCIIAT( JWAVE + KLAST »yNUMBy PRINT) 
DIMENSION YO(400) 9RO(400) 9SC072) 
COMMON /NO4/Z(400) yAAMP (400911) »BAMF (400911) 
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8 CONTINUE 
TOTAL=TOT/SUM 
WRITE (6210) TOTAL 
10 FORMAT(’ “’»9’FOR FIRST 6 WAVES VARIANCE EXPLAINED’ 


LrF7,.3) 

RETURN 

ENT! 
C 
CXXKKKK FORECAST TIME WAVES IN AAMFLy BAMPL» 
2 AND DIFFERENCES  XxXXKXKX 


SUBROUTINE FCSTCNWAVE »ROrL e TEL » JT »NUMBy VALUE » 

LAMFPLOsFRINT) 

DIMENSTON VALUE (400) sACIS1) 

COMMON FREQC(1S1) sAMPL (C151) sPHASECIS1) xAACISI) 
DATA FI/3.141593/ 

Q=2.KXPT/3460.4 

X=2.KPTRTEL 

KLAST=NUMEFL 

K=0 

C *x BY-FASS THIS LOOF IF CLIMATOLOGY FORECAST WANTED 

TF CABS CARS(CPRINT)-S3.5)-LT.0.1)960 TO 8 

NO 7 T=isNWAVE 

ACT =2.xXAMPL CT) 


K=K+ i 
C x CHECK TO SEE IF THE AMPLITUDE IS LARGE ENOUGH 
C TO BE USE ? IF AMFL LESS THAN AMFLO THEN NOT USED 
IF CAMPL(I).-LT«AMFLO)GO YO 8 
7 CONTINUE 


8 UWRITEC6r10)K 
10 FORMAT(’ ’»’ THE FIRST ’*IS*’WAVES WILL BE USED’) 
JORG=1+ (NUME/2) | 
C *k JORG=ORIGIN 
nO 3 JJ eKLAST 
SUM=RO 
C xx CHECK FOR CLIMAT FORECAST 
IF CABS CABS (FRINT)-5-5)+LT.0.1)9G0 TO 6 
JN=J-JORG 
T=FLOAT CIN) 
C ok ALL FORECAST WAVES TO MEAN OF AMPLITUDE 
n0 4 I=iyk 
SUM=SUMEA CL) XCOS (XKXFREQ (I) KT+QXFHASE(T)) 
CONTINUE 
VALUE (J) =SUM 
CONTINUE 
RETURN 
END! 
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C 

CKKKKKK FORECAST REPORTS  X#x4KX 

C 
SUBROUTINE FCIIAT<( SUAVE» KLAST yNUMB ys PRINT) 
DIMENSION YO(400) ¥RO(6400) S072) 
COMMON /NO4/Z(400) yAAMF (400911) 9 BAMF (400911) 
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8 CONTINUE 
TOTAL=TOT/SUM 
WRITE (6210) TOTAL 
10 FORMAT(C’ ’»’FOR FIRST 6 WAVES VARIANCE EXPLAINED’ 


LrF7.3) 

RETURN 

ENT 
C 
CXXKKEX FORECAST TIME WAVES IN AAMFLy BAMFL» 
: AND DIFFERENCES  XxXXKXKX 


SUBROUTINE FCST(NWAVEyROrL DEL » JI »NUMB» VALUE » 
LAMELOyF RINT) 
DIMENSION VALUE(400) sACI51) 
COMMON FREQ(154) 2AMPL (154) sPHASE CISL) rAAC1S1) 
HATA FI/3.141593/ 
Q=2. KF 1/3460. 
X=2. KPT XDEL 
KLAST=NUMEEL, 
K=0 
C *k BY-PASS THIS LOOF IF CLIMATOLOGY FORECAST WANTED 
IF (ARS (ABS (FRINT)-5.5).LT.0.12G0 TO 8 
10 7 LT=lyNWAVE 
ACI) =2.kAMFL (1) 


K=K+ 4. 
C %*x* CHECK TO SEE IF THE AMFLITUDE IS LARGE ENOUGH 
C TO BE USE : IF AMPL LESS THAN AMFLO THEN NOT USED 
IF CAMPL(I)-LT«-AMFLO)GO TO 8 
7 CONTINUE 


8 WRITECG210)K 
10 FORMAT(’ ’*’ THE FIRST ‘’»ISr’WAVES WILL BE USED’) 
JORG=1+ (NUME/2) ‘ 
C xk JORG=ORIGIN 
NO 3 JeJIyKLAST 
SUM=RO 
C *’k CHECK FOR CLINAT FORECAST 
IF CABS (CABS CFPRINT)-S.5)+LT.0.1)9GO0 TO 6 
JN=J-JORG 
T=FLOAT CIN) 
C *k AD FORECAST WAVES TO MEAN OF AMPLITUDE 
00 4 IT=isk 
SUM=SUM+A CT) XCOS (XXFREQC TI) XT+QXPHASECT)) 
CONTINUE 
VALUE (J) =SUM 
CONTINUE 
RE TURN 
ENT! 


Wa D 


C 

OOKOOK FORECAST REPORTS  okx0KxKx 

C 
SUBROUTINE FCIIAT( JWAVE » KLAST »yNUMB ys PRINT) 
DIMENSION YO(400) »RO(400) S072) 
COMMON /NO4/Z(400) vAAMF (400911) » RAMP (400911) 
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DATA TORT /6«2831853/ »Fi/1.388S889E-02/7 
C xk FI=FUNDEMENTAL FREQUENCY 

READCO»sLIKLASIs MA 

1 FORMAT (215) 
IF CKLAST-EQ.KLAS1)GO TO 2 
WRITE CS 3) 

3 FORMAT(’ “s’TREND ANT KLAST [10 NOT MATCH’ ) 
STOF 


REAIICS »4)> CVYOC LD)» T=1yKLAST) 
FORMAT CS rTP Ay 3) 
010 S T=NUMBsKLAST 
C xx ADD DIFFERENCES TO TREND TO GET FORECAST 
C MEAN OF THE REFORT 
ROCTIRYOCTIOFZ¢01) 
CONTINUE 
IFRINT=IFIXCFPRINT) 
WRITE C1391) KLASTs IPRINT 
NO 6 IT=NUMBsKLAST 
C xk THIS WILL FRINT OUT FORECAST MEANS OF THE REFORT 
C IF /FRINT/ GREATER THAN 7.5 
LFACARS< PRINT) St se oO GURILO RTS 
WRITE CG»10)ROCT) 
10 FORMAT(’ “’9’ROCI)=% 9 F7.3) 
12 CONTINUE 
C **k USE FORECASTS OF AAMFL AND BAMFL TO 
C FORECAST REFORTS 
NO 7 Je1ls72 
SUM=ROCT) 
X=FLOAT(J-37> XFL XTOPT 
NO 8 M=1»JWAVE 
D=FLOAT CM) XX 
C=COS(T) 
ST=SIN(D) 
SUM=SUM+2 + KAAMF (I 9M) XC4+2.XRAMP CT yM) xST 
8 CONTINUE 
§(J)=SUM 
7 CONTINUE 
C *xk STORE FORECAST REFORT 
WRITE (1329) (SCJ) 9J=1772) 
9 FORMAT(’ %93¢/24F4-0)>) 


C *x GET TREND FROM STORAGE 
a 
4 


ay 


& CONTINUE 
RETURN 
ENT! 
C 
CKORKKHOKK SMOOTHING SFECTRUM XXX KKKX 
C (SEE CHAFTER 3 >? 


SUBROUTINE SFCSM(JJ,SUM»sFRINT) 

DIMENSION SMFHC1S1) sSMAACLTS1) »SMFRCLS1) sACI SL) 
DIMENSION OLDUR(151)*ASQ(151)*BSQ(°151) »AABSACIS1) 
DIMENSION B(151)sBRARSECIS1) 

COMMON FREQ(151)2R(151) sQC151) rAACIS1) 
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DATA TOFI/4.2831853/ 
T=TOFI/340,. 

JH=JJ-1 

JP=IJI-2 

FREFARE DATA FOR USE IN SMOOTHING OF TIONS 
NO 1 J=1sJH 
FHASE=Q( J) XT 
ACJO=RCJ)XCOS (PHASE ) 

BC J)=-RCJOKSINCPHASE ) 
ASQACIJORACII KA CI) 
RSQ(J=BCJOKERCY) 
AABSACJI=ACI) XKABS CAC I) ) 
BARSBCJO=BRCJKXABS CRC) ) 
CONTINUE 


CHOOSE OF TION 

IF CABS CABS (CPRINT)-3-.0),.LT.0-1)G0 TO 300 

IF CABS CABS CFPRINT)-2.0).L7T-.0-.1)G60 TO 200 

IF CABS CABRSCPRINT)—-1.0).L7.0.1)G0 TO 100 
WRITE (G92) 

FORMATC’ “9 °XXXNO SMOOTHING OF TION USED Xxx’) 
RETURN 


WRITE (49301) 

FORMAT(C’ “9% SMOOTHING OFTION NUMBER 3 USED’ ) 
AR=(AABRSACLIFAABSAC2))/2. 

BR=(BARSR( LI +RABSER(2))/72. 

SMFPHC1)=ATAN2 (-BRe AR) ZT 

SMAACL)=ABRS CAR) +ABS CBR) 
SMFRCL)=CFREQCL) XAACLITFREQ C2) XAAC2) I /CAACLITAAC2) 
NO 302 J=2y JF 
AR=(AARSA(J-1)+AABSAC JI +AABSACJF1)) 73, 
BR=(BARSB(J-1)+BABSR(JO+BABSBCJt1) 73, 

SMPH CJ) =ATAN2 (-BRy AR) /T 
SMAA(J)=(ABS (CAR) +ABS (BR) )X2.73. 
SMFR(J)=(FREQ(J-1) XAACI- 1) +FREQCI) KAAC I) FFREQCIt1 DX 
LAACIFL) YS CAACI-1IFAACIJIFAAC ITI) 

CONTINUE 

SMAA (JH) =AA (JH) 

SMFR (JH) =FREQ (JH) 

SMPHCJH)=QC0JH) 

GO TO S00 


WRITE (69201) 

FORMATC’ ’»’ SMOOTHING OFTION NUMBER 2 USED’ > 
AR=(ACLI+A02))72. 

BR=(BCLI+R C02) )72. 

SMPHC1)=ATAN2(~BRyAR)/T 
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100 
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300 
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SMAA(C1)=AAC1)+AAC2) 

SMFRC1L)=(FREQ(1) *AAC1)+FREQ(2) KAAC2) )/SMAAC1) 
SMAAC1)=SMAA(1)72, 

BO 202 J=2» JF 

AR=CACI-L) FAC J) FAC I419973. 

BROCE Cia bR +E Cee) 7/3 

SMPHCJ)=ATAN2 (-BRe AR) /T 
SMAACJIB=AACI-LIFAAC I) +AAC IHL) 
SMFRCJD=CFREQCI-1) XAACI-L)+FREQ CI) KAAC JO FFREQCIt1)% 


ZAR CIFIID/SMAAC I) 


SMAACIJO=SMAACIJIO/3. 
CONTINUE 
SMAACIJH)=AAC JH) 
SMFR( JH) =FREQ( JH) 
SMPHC JH) =Q (JH) 

GO TO 300 


WRITE (S101) 

FORMAT(’ “’s’% SMOOTHING OFTION NUMBER 1 USED’) 
AR=CACLIFAC2QII 72. 
BR=CRCL)F+RC2))72, 

SMAACL)I=CARKARTERKXBR) K2¢ 
SMFR(C1)=(FREQ(1)XAACLITFREQ(2)XAAC2))/CAACLITAAC2) 
SMFHCLI=ATAN2(-BR se ARD/T 

10 102 J=27 IF 

AR=CACIA-LIFACIJOFACIFIO IZ Se 
BR=(BCJ-L) FRC JBC S12 0/35. 

SMPH( JO =ATAN2 (-EBR YARD ZT 

SMAA (J) =CARKXAR+BERXBR) 2, 
SMFR(J)=(FREQ(J-1) XGA C I~ LI FFREQ CJ) KAAC J) FFREQCUt1)X 


SAAC IFLIII/ZCAACI-“L) FAACUD FAA CIEL) D 


CONTINUE 
SMAAC JH) =AACJH) 
SMAAC JH) =AAC IH) 
SMFR (JH) =FREQ( JH) 
SMFH (JH) =QCJSH) 


CONTINUE 

SRMAX=0.00 

MOVE SMOOTHED IATA TO COMMON BLOCK ELEMENTS 
NO S01 T=1+JH 

AACT=SMAACT) 

OLDR(I)=SQRT(SMAACT)/2.) 

SRMAX=SRMAXTAA CT) 

Q¢L)=SMPHCT) 

FREQ(I)=SMFRCI) 

CONTINUE 

WRITE (62502) 

FORMATC’ “710X2 ‘FREQ’ 910X»9 “AMP ’ 912X% “OLD! KR’ 2 8Xy 
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IF (SRMAX.LT+0.0001)SRMAX=0.0001 

SMOOTHING FINISHED--- NO VARIANCE ADJUST 
00 S03 T=1»JH 

AA CT )=SUMXAACT )/SRMAX 

RCT=SQRTCAACT) 7/2. 

IF(I.GT.40)GO TO 503 

WRITE (S504) FREQ( I) es RCI) sOLDR(I) sQ¢1) sAACT) 
POURM OWS oe Sk baat? 

CONTINUE 

WRITE (49505) SRMAX 

FORMAT(’ “s “TOTAL SMOOTH AMPLITUDE SQUAREL 
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APPEND IA st 


PROGRAM PARAMETERS 


IWAVE = number of strong waves that are to be listed. 
These waves have been sorted according to the 
magnitude of the amplitude. 

PRINT = parameter that determines smoothing option to 


bes useds (see. bellows descr piaonye 
L = number of forecasts to be made. 


DIV(=JWAVE) = number of space-waves used to forecast 
heights. 


NWAVE = number of time-waves used to forecast amplitudes 
for space-waves. 


NUM = NUMB = number of elements in the diagnostic data 
SHOE Sc 
wy 
RMS = (Variance) “<"5 The Variance’ of differences from the 


trend 


Hie-wstanuing DOunt sOfedlacnos |) Caddaamseis 


Values for PRINT Parameter 


If print is negative, then the unsmoothed spectrum 


JESMEILOLL MAAS CEG, 


Absolute Value of PRINT 


[0 252 no GSMOOt NIN GROUtN ON MUS CGmD Uisms Ob CUms Enon 
waves not printed out. 
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20m Onl se smocthingmont ron: Noe 2 susede 


GT 3.5: no smoothing option chosen, therefore, bypass 
smoothing subroutine SPCSM 


= 4.0: persistence forecast made (see verification 
program listing in Appendix I) 


5 OD climatology forecast made 
note space-wave pattern is constant but 
means still adjusted for trend 

= 6.0 value that can be used when verification 


CoGirekn BSe) Meee. 


GE Os IDrint out forecas tumeans 


GT 9.5: print out estimated AAMPL for forecast period 


LT j= less than 
GT = greater than 
GE = greater than or equal to: 
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APPENDIX J 


GISTING OF VERIFICATION PROGRAM 


Files Needed for Input 


Number 
1S —“ROrecadast MencgntenMeDOrLs TroMmeramiec Of the 
program. 
5 -Actual height reports. 


Note: Persistence forecasts are generated by lines 
42 to 47. 
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CG 
C XXK VERIFICATION FROGRAN x*x*x 
COCO GCC CORIO OOO OR OK KOK ak 
C 
Bs 
DIMENSION SC72) s¥(72) e072) sA072) sR (60) sRMSEE(40) 
INTEGER HOV( 40972) »yFHOVC 60972) 
READCLS» 1 KLAST» LERINT 
1 FORMATC2I5) 
TOTAL=28.0 
VO 2 [=1%60 
NOZO=0 
Loa) 
TIME=FLOATC(CIO)/2, 
WRITE C6*78)T0*TIME 
78 PORMAIG 604s “FOUREGAST’ylSs"=’s5Xs*T0t7 sbS. is” DAYS") 
REACT Ss 9) CSC ee da1s72> 
G) BER MA TG Se3t 72454 .0)) 
READS S8d CY CJ 9 J=1972) 
8 FORMATC2°024F3.0°8X/) »¥24F3.0) 
XX SET UF HOVMOLLER DIAGRAMS 
HOV= ACTUAL HEIGHTS 
FHOV=FORECAST HEIGHTS 
Wad pie TB cos ree 
HOVCTy»yLO=IFIXCYCL)) 
FHOVC Pei etr iA eS th 
IFCY(L)Y.EQ.0.00)NOZO=NOZ041 
12 Cea Or Aha E TE 8G an he 
ACLI=YCL) 
Lo GONTINGE 
C ** CHECK FOR MISSING DATA 
LECNGZ OSE Tae COV TOI LG 
WRITE (S6917)N0Z20 
17 FORMATC’ ’9% DATA MISSING NO ZEROS=% ¥IS) 
RiGlw=o9 94 27 
RMSEECI)=-9.99 
TOTAL=TOTAL~-1. 
GO TO 2 


oa} (eat (ep 


cma] 1 Lge 


**k CHECK TO SEE IF FPERSISTANCE FORECAST WANTED 
16 IFCIFRINT.NE.42G0 TO 21 
NO 21 Jelr72 
S¢( J =ACI) 
FHOV( Ty JO=IFIX(SCI)) 
21 CONTINUE 
C **k PRINT HEIGHTS 
WRITE CG919) 
19 FORMATC’ ’+s*’ FORECAST REFORT = S(lo* > 
WRITE (6977) (SCJ) re J=1972) 
77 FORMAT(’ ’93¢C24F4.0/791X)) 


L Gia 


we Hor pokes 
(ove dont: * Spaibety 
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WRITE (C618) 
FORMAT(’ “s’ ACTUAL REPORT = YC(L)7) 
WEEE CGriee ) CY Cl) s led to) 


C X#kKK CALCULATE STATISTICS FOR THIS FORECAST 


C 


im 6 


13 


r 


a 


AND ACTUAL HEIGHT DATA AT THIS TIME 
COF=99, 
50=0.00 

SMAG=0,00 
YMAG=0,.00 
SUOTY=0.00 

SUM2=0,00 

SUML=0.00 
SUMN=0.00 
NO 3 K=ly72 
QCRI=SCKI-Y CK) 
SUM=SUM+A8S CICK) 
SUML=SUMIF0N CK) ATICK) 
SUM2=SUM24ICK ) 
SO0=S CK) +50 
SUR ta Ok CR pe STO 
YMAG=YCKIRXYCK) + YMAG 
SMAG=SCKIXSCK) + SMAG 
CONTINUE 


SMAM=SUM/72. 

DUMEAN=SUM2/72. 

SMSE=SUMI/72. 

SMYMN=YMAGXSMAG 

Tecate. 0.0001) GO TOs 
COF=SrOTY/SQRTCSMYM) 

SO0=S80/72. 

WRITE (4912) COF 9 SO 

FORMAT C’ 49 “COF=%¥F7.395X» “SO=% 9 F723) 
WRITE (4694) SMADy MEAN» SMSE 

FORMAT(’ “» “% MEAN ARSOLUTE UIFFERENCE=% »F7.3%73Xy» 


1/MEAN DITFFERENCE=’ »F7.393X»’MEAN SQUARE ERROR 7 


y “ESTIMATE=% sF 7.3) 

RMSEE (CI)=SMSE 

YO=SO0-DMEAN 

SY=SQRT CYMAG/72.-YOXKYO) 

SS=SQRT (SMAG/72.-SOxSO) 
RCOF=(S0HOTY-72.*SOKYO)/072.¥SYXSS) 
R¢T)=RCOF 

WRITE (6913) YOrSY*SS»RCOF 


FORMAT(C’ °s 7YO=/ sF7. a9 5Xe/SY=" sF7. 32 GXs “SS=" eF/e Sy 
SoxXs R=’ sF 7.3) 


CONTINUE 
(CONTINUED ON NEXT FAGE) 
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105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
hi? 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
134 
132 
133 
134 
135 
136 
E37 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
ie. 


153 


C 


C XXXK FRINT OUT FIRST 28 FORECAST VALUES 


C 


C 
C 
C 
C 
C 


23 


[5 
aoa 


FOR KyMSEEsAVERAGE Rs ANT AVERAGE MSEE 
WRITE (46924) 
SUM=0.90 
SUM2=0.00 
TO=0 
WRITE C6923) TQrR¢C1) »RMSEE C1) 
NO 22 T=2929 
TO=T-1 
WRITE COs 233 TO eR CID eRMNSEEC TD 
FORMATC’ “9TS9F?7.39F 7.1) 
TPCARS (RCT )-99,.99)..7T.1,.00I60 TO 22 
SUM=SUM+R (CT) 
SUM2=SUM2TRMSEE (CT) 
CONTINUE 
FORMAT (417 »3X9°I% sSX2’R% s SX» “MSEE ” ) 
SUM=SUM/ TOTAL 
SUM2=SUMN2/TOTAL 
WRITE (S925) SUM eSUM2 
FORMATC’ “’9 “MEANS ’sF7.39F 7.1) 


*KKK PRINT OUT HOVMOLLER DIAGRAMS **x 


31 


ACTUAL DATA 

IF (CIFPRINTXIFRINT)-LT.9)GO TO 30 
WRITE (4928) 

FORMAT(C’1’9 “FART 17) 

NO 26 T=1y2t 

WRITE (4927) CHOV( Ts J)» J=19346) 


7 FORMATC’ %93613) 


CONTINUE 
WRITE C6929) 
FORMAT(’1’9’F ART 2%) 
LO 30 T=1%28 
WRITE (6927) CHOVC I» J) 9 J=37 2772) 
CONTINUE 
FORECAST DATA 
WRITE (6928) 
VO 31 T=1s28 
WRITE (6927) (FHOV( Ly J) yJ=1%36) 
CONTINUE 
WRITE C&¥29) 
VO 32 T=128 
WRITE (6227) (FHOVC T+ J) 9 J=37972) 
CONTINUE 
STOF 
ENT 
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TEST DATA PROGRAM 
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CKXKKK FROGRAM TO GENERATE TEST DATA *xKxXxX 


C 

C 

C 

C 

C 

C 

C 

C xXx 

C XxX 

C xx 
13 

C 

C 

C xx 

C 

C 

C xx 

C XxX 

C xx 

C xx 


S=GENERATED HEIGHT REFORT 

A=COSINE COMPONENT OF WAVE AMPLITUDE 
R=SINE COMPONENT OF WAVE AMPLITUDE 
FyC=AMPLITUDE/2 OF THE SFACE WAVE 
QyCQ=THE CORRESFONDIING FHASE 


DIMENSION ROC400)*S¢72) ACL1) es BCLL) eC C11) sQC1d)» 


IF C11) »CQC11) 


INTEGER 7 

VATA. URI 7 G5 2as3.167 

F=TOFI/730. 

TREND) COEFFICIENTS INATIALIZED 
ROM=51.0 

G=-4,0 

H=~-11.0 

Z=TOFI/360,. 

NUMB=90 

THE CENGTHOOF THE CONGEST: Wave 
CUMB=FLOATCNUME)-460, 

UTS USED TO UETERNINE THE ORIGIN 
U=CUMB/2.+1. 

WRITE (7s13)NUMB 

WRITE C8s13>NUMB 

WRITE C9s13)NUME 

FORMATCIS) 


GENERATE THE TREND 

tO 1 T=1»NUME 
ROCT)=ROM+GXCOS CP XT) + HKSINCRXT) 
X=(T-U) /CUME 


GENERATE SPACE WAVES 

DO 2 M=1y114 

V=FLOAT (CM) 

SET MEAN VALUES FOR ACM)» RCM) 
ROMB=1.5/¥ 

ROMA=-0.8/V 

SUM2=ROME 

SUM1=ROMA 

All TIME WAVES TO SFACE WAVE AMFLITULIES 
10 3 I=1r11 

W=FLOAT (I) 

SET SFACE WAVES=0.00 
C(1)=0.00 

F(1)=0.00 

Q(1)=0,.00 
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KK 


hh 
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CQ¢T)=0,00 

FICK WHICH WAVES WILL BE NON-ZERO 
TPCl-ER M)GO TOS 
TROIS EG. ekrnooGO Tors 
IFCI.EQ.(3*M) GO TO 5S 

GO TO 4 

SET THE VALUE OF THE NON-ZERO COMFONENTS 
FCT)=0,9/W 

QCT)=40, *XWKZ 

CQCTI)=-10. kWKZ 

C(I)=0.7-0,03kV 

CONTINUE 

E=TOPIXXXW 

CALCULATE THE SFACE WAVE COMFONENTS ACM) »B(M) 
SUM2Z=SUM2+2.XF CI) XCOSCE+Q¢(1) > 
SUM1=SUMI+2.*C (CT) XCOSCE4+CQ¢(1)) 
CONTINUE 

ACH) =SUM1 

BRCM) =SUM2 

CONTINUE 

WRITE (7212) (ACJ) + J=1911) 

WRITEC 8912) (RC J) sJ=1911) 

FORMAT C11F 7.2) 


COMBINE TREND ANDI SFACE WAVES TO GET REFORT 
DO @ ea le= Liv 7 2 

SUM3=RO¢T) 

tti=TOPIXFLOAT(L-37)/72, 

0 6 M=1lyii 

DiM=DKFLOAT CM) 
SUMS=SUMNSt2.¥*ACM)XCOS CIM) +2.XB0M) XSINCIM) 
CONTINUE 

S(L9=SUMNS 

CONTINUE 

STORE THE REPORT 

WRITE CLOs1 5) (°S¢(L) »L=1y24)9T oT 

WRITE (10916) ¢S(L) »L=25%48) »T 

WRITE CLOxyI 7 (SCL) yL=499 72) yROCT) 9 X 
FORMAT (24F 3.07214) 

FORMAT (24F 3.0914) 

FORMAT (24F3.092F 4.1) 


CONTINUE 

STORE THE TREND 

WRITE C9914) (ROCT) » T=1 »NUMB) 
FORMAT Gs Libv. 3s) 

STOF 

ENT 
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Table Li: Results using other various, Cech- 


niques for smoothing. (Used on test sample) 
Method: 

1. Smoother Option No. 2 with variance adjustment. 

2. Smoother Option No. 2 with no variance adjust- 
ment 

3. Smoother Option No. 2 with no variance adjust- 
ment and the smooth mean square of the time- 
wave withewaveonumber WwW ecalculated@=py (0 + >AA()) 
FORA Gai) if Se 

4. Samecas method 2, except the frequency is not 
smoothed (i.e. not weighted) 

5. Same as method 1, except the time-wave with wave 
number 1 1s not smoothed. 

6. Same as method 4 except the time-wave with wave 


number 1 is not smoothed. 


The value vr in «thls ADDENCine 1 SeliecOnne- 


lation coefficient (i.e. r = R) 
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(Southern Data Set) 
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Results Using the Southern Data Set NUMB = 90 
and JWAVE = 8. 


Wectlavich syle & 


Results Using the Southern Data Set NUMB 240 


it 


Thevalue ur i) CHES BADDeNG lees ts. LNe SCOT nes 
lation coefficients (ieee =~ ts) 


NOTE: Forecast numbers are listed to the left of these 
tables. Forecast number = 0 as forethe last diagnostic 
report, and 1S 0b considered siecnesMeanmat= thes bottom 
Otemther tables 


vat “ 


; 7 
1) anf Ne 


Cee, wnat wie 


@ 
Aruve Bre gabe aaiy a Ran 
: dl 5” ch ia Wis he 


: : 43% r a 
% = 
7 ; 
—_ : 
oa 
af 7 4 2 i) 4 ae 
, ‘= A el 
Nd ' 5 _ 7 ‘ 
agape, UMS, aT PA eet bohen TT) eae 
: rt, 4 oy Boll Pe ee F519" otra _ 
as . 


wi bi pe gil pra. Fredy! pediath ay 

cor 4a i P Meet egser0%) vt 
Aan ebt at tava nN) A had a 
9 . . fy _ oy Wey 


153 


Lie ov LU 
UOT Se 4 "0 
ANE ONE 8 
O ieee 9.980 
(ee Blake) 18 
de ee 18 
G6 ede 18) 
ay LoL. 0 
c8 eZ 0 
BL oye 18) 
ae 98020 
£9 2s) 
as sue ae 
09 Gla 0 
Users | oaeU 
86 965 0 
VO eo G20 
Lee. Ga 0 
Colma Se (0 
God oo Va0 
NS) Gos A 
7 pee Os 
SO Tee 7 0 
Vo bees 7 20 
86T ver'0 
Gola 65.0 
JISW A 
Z=J4AVMN 
cae ON 
AIYUJOOWS 


HU Leo © 0 
Oi [auc 0) 
CS tee LU 
ls eves ie 
S61 v6c 0 
Uo lee Goes 
LONG (oat aD) 
CME ABU EE 18 
HO Tors eee 
CU eco a0 
SOME Ase 1 
€6 685 0 
ob ESB) 
08 009 °0 
CU ee G0 
ql AR NSGH a0) 
Cee 
AG Re 48) 
Se Basis) 
16 (AGL 38 
vl ULL 0 
ie 8280 
vs zisus. 1) 
aG rates: 10 
8S Coe 
6€ eco U 
4ASW J 
c=3AVMN 
YYOOWS Uf} 


ede 


0 


Ce Se (8) 
OL Oy. 10 
fis ouky, 0 
06 Sree. 1) 
26 Cie) 
G8 LoL 
tl Osea) 
vs Stowe 10 
19 Cie 30) 
she Shjays, 8) 
CG Cea 
a cio 0 
Oe chairs 1) 
OV cvs 0 
GL 6090 
wa Sis 0) 
6L s/s) 
€6 sie 0 
L6 eae 
HOM oe 
96 6£9°0 
€6 002 °0 
DORbeaeee OU 2en() 
Sah Maite (0) 
Vice oon 0 
Aah eho) 1G 
JASW dl 
T=3AVMN 

c¢ ON 
4ayuzOOUWS 

8 = JAVMC 

6 = GWNN 


98 €09°0 
v8 Siege) 0) 
ExG Aree 8) 
WA IE | TRIS AG) 
Go Oa 
Ee yale 8) 
POS Suse de 
ae Growme 
08 €89°0 
08 139-0 
98 eohs) 18) 
EL 169° 0 
v9 pOZ°0 
L9 CUO 
68 hye) 
He Higiey 18 
CO Tec ood) 
WALA — [efeiSe 49) 
Sa Gee 
66 WISE 
98 9GL0 
6L ct eea) 
6L ahi 8 
88 p8l°0 
88 964.0 
rare 998 0 
JASW J 
L=3AVMN 
YYCOWS Uf} 


96 

ae! 
c8 
Sut 
Gus 
90T 
SOl 
88 
v6 

96 

pol 
fea yl 
vel 
Saleh 
Ol 
Shine 
On 
Skt 
Sol 
06 
BL 
6% 
Tv 
Ge 
ue 
0 


44SW 


BIUaYSLSUaq 


Eisen (0) 
758 0 
9€8°0 
CLiwau 
Cas) 
OL LO 
Olan 
eres 10) 
Gib} 8 
964.0 
Chs Aaa) 
sv L40 
psy Ae 10) 
Ca) 
sto 7ee 18 
Shine 48 
BG L400 
Loe) 
One Lea 
S080 
0€8°0 
pL8 0 
Eber ie) 
676 0 
EG6 0 
000°T 

J 


is 
St 
8S 
vl 
c8 
oe 
69 
90 
Es 
oS 
oa 
BL 
6€ 
LE 
60 
0S 
v9 
(ae 
Ls 
08 
as 
ie 
Vi 
62 
L8 
Gil 


JASW 


6180 Ge 
6€8°0 vc 
G6 oan) EC 
Via cc 
ence 0) We 
itl eo ( 8) O02 
Ceae 6) 6l 
928 °0 oa 
608°0 1a 
S13 30 oF 
028 °0 GT 
erent vl 
8Ps8 oO El 
028 °0 ol 
sibrss 18 Tt 
323 0 OT 
LSLGU 6 
Corey) 8 
Oo LO Me 
Sue 18. 2 
alae Aw) S 
ese v 
v8 °0 € 
0€8°0 é 
028 °0 if 
£06 °0 0 


A 


Aboyoyewr ly 


19S e2eg UvaYy INOS BHuLsN sz[Nsay 


“TOW PLGeL 


ic? i oo Cs) 
met a' 


“i 
®@ 


weiss 


= 


—_ 
yo =O 4 


eo we 411 Fe ! 
2 er. 
smh & oe 2 
ihe Wee > 
9 69 4 Gp 44 10) & 


rm 


- 


eo on. woe 
ef sfed sovSitucet betel prfucee : 
: ~~ : ea 7 
*cahpest”— <—25ep7de+94 
Puy TARN : 


s 


ier 45. 1 Fe 


= 
of 


(a Aam 


i i ie) 


=—— & 
& S&S 
rr 
ots 
/ F 
ea 
end 7 


eca3o— a pnd .1 at) = 
— @F% 26 es a té 
—— 28.2 <3 tse .¢ . 
5 ie~ s)¢e¢m ; 68 
28 NV cl ce el=c9 ry pode 
£336 Bt = o1eet ¢ 7cv.8 | 
& a GY reoul DS se % 
tbe 4 Fo) Ott7.0 \c. ; atty, 
t (e208 ist. wet S- ‘62% 
$—- ff? .g Ctl: fet. ba. 73% 
cau ott ee a < 
5-~e2.6 Gat itt.9 os 
> ie 1 ett -stt.2 te 
a4 e2*-_0 >t a 
a =f 
, Dp D> 2. oe 7 
L842 .= of ;.§ 3 
fe »s if. te 
aa m f 256 .f sf 
z fe st ar) 
é 7 * oll 2% 
; ec a ics ..0 98 
7 wy Oi 7: ag 
ie ,< ae ‘e 
: » a a 2 $3 


- 


154 


Gel 


SHEA 
Gel 
80T 


sas eh ad Osa 
oe) a 
c0L°0 8ST 
é92£°0 OST 


O€S 0 


686 0 
Boca) 
O€€ 0 


SL 


iat 
v8 
OL 


Ci tan c6 
hye 18) vol 
Eerie 10 vol 
Lo Gad iS) 


sine) 1) 


CigceaU) 
tse, (0) 
GUG <0 


G6 

est 
aa 
vol 


,3u09 


€08°0 
67L°0 
ed 220 
908 °0 


ie a Suck 


v9 86L°0 
98 Daan 
08 0S ae 
eg Bus <0 
"TW 9LGel 


NVAW 


8c 
Le 
9¢ 


* "sna ay oe 


i te 
SN. “gant £5 ae tre est.0 


23%. a - 7 7 a2e 0 29 ce .0 


S bd 7 
7 e 
7 7 : 7 
4 7 — 
' ; = “ - - 
7 _ 
¥ : : = 
7 7 2 - - = < 
ee =~ “a we 
ew, ~ =e od 
= £ 7 y 7 
mae. 
= ; 
= 
7 : kr _ 
7 7 : r » 
ia P [ - 
a * * 
: 5 sh 
x 2 ‘ be x “ - -_ 
_ -_ 


sy 


vais ASRS 10 
JEM ASHEN 
OO ee bo ert) 
Go Le ve 0 
Geer le ce.U 
vee Bear 0 
O42¢e o09¢ 0 
Sa 103 a8) 
rele 777.0 
9c Ga 0 
SOL 1os 0 
Vetere oo 80 
66 63s 0 
€8 C2940 
68 £09° 0 
Cee cca 
AVES 30) 
wodh ~ ase 
[oth ee) 
Cole 2 Ly..0 
Sy eas 1) 
Pa a0 
S6 BOLO 
9L ned AW 
cg 7953 °0 
ge O Lea. 
ve SG s 0 
4JISW J 

G = JAVMN 


pay.yoowsun 


Ssvleeg9 0 
Hse “SCRE Ta) 
Ole 61 50 
Cra le JON 10) 
Oodle eos 0 
EL igs7S9 30 
Lo 9 O0'0.0 
9/ vl 0 
(Go SV 
(Re Oe 18 
Cc om 
Power 6 Oa 
8b 0280 
Bam Cee 18) 
We) fase) 19 
ete coven) 
Cre eo ae) 
Cel (ein 0 
Gul suv 0 
(sly aetsaie 16 
Cia Week) 8 
(eS gas 18 
AC Diy Jae 
eee Seah 16) 
Av AON 19) 
VS ISO LGs ty) 
Cee Sa) 


JISW A 
v=3AVMN 

¢ “ON 
AdYPOOWS 


qued 


se tease 10) 
WELT Seca 0 
Gol. tres0 
Ost 9 coe0 
Glgeec Go.) 
Gao cecu0 
UR Jeera 8) 
Oa gee eo 
Vitae Le 0 
66 7 GeO 
56 See 0 
0 Cena ce 
06 Oe) 
ae v9 0 
08 ess) 0) 
Se en) 
Zi bes Se 
Pe 729 20 
Oc led cae 
oSt O27 0 
eee AOS Se 8) 
aha“ ON8) Ae) 
26 OLZ 0 
ae 882°0 
eo sre: 1 
85 swlten 10) 
LE 676 0 
4JASW J 

fb = JAVMN 
pay zZoOOWws UY 
= JAVMC 


GWM N 


9p1 985°0 
8Z1 169°0 
PEL 929°0 
gzt rs‘o 
es Se 
801 $69°0 
ES te 
€9 €08°0 
Z2y = 080 
6p 0£8°0 
ede AIA 
re) Se 
Ce Ee 
99 §99°0 
6£ ¥09°0 
OZT 62r°0 
Zo Sew 
prt 66€°0 
ZGtelO7 0 
eGi -61b 0 
LGieelc uo 
9pl $90 
Lvl 9L0°0 
ae 
vik evn 
y8t o8sr'o 
8ST 019°0 
FERS. J 
€=JAVMN 
€ ‘ON 
AIDYPZOOWS 


19S eeqg uuauynos BHBulispn szLNsay 


6G e yoor0 
CcCheeer Ou G12) 
Cie Lec vy o0 
Gh 6 690 
CCC oO 
Oy eee ig 
96 COLO 
08 622-0 
vS (ese 10) 
8g aMEyE 19) 
8S G9L°0 
GS GUL 
60 Olecan0 
vv ye 0 
LS Oye ye 
URGibe Sele 8) 
LOE USse0 
ev orc aa) 
Seah © eS a0) 
Le eee 
ce eo Ca) 
Oe 20 
92 ~ 19S 9 
GMa me ASSIS 16) 
ony esa secr 148) 
Tee Oe GeO 
EG 2bo"0 
J4SW A 
CAAVMN 

c °ON 
ABUPOOWS 


B91 eco Came C 
Oot Soe s0aee ac 
OV oc Oe eer 
OG a vce ec 
Got (OLcs0 sec 
C0 leeeL Ue Ome 
LISS Can 0 ae 
PG Laue S0Va weet 
SHOE eariiecy a0) thi 
66 IESE AOS eel: 
66 ese at) 1S) 
COL 7 Sa eee 
66 8Z9°0 wl 
o8 909°0 €T 
68 v6g°O el 
Oe ss el Uk 
CC Pees me 
ORO ei Sp 8) i) 
ics ee Cea 
6bl vp6r'o L 
SAI calonsy: eh 
96 Ob9e0 ms 
SL STE AVE: 51) 
od ay th os 
vS Sashes (bh 4 
Sls 988°0 | 
site 2S Oa OmaeO 
4JISW J 
€=4AAVMN 
poyzyOOWs Uf} 
"TW 9LGPL 


Kg j 


oe 
- 

- _ 
2 ; - _ 
“a -@) “a= = om SIs 
Awe lte« et «be = yer . oD he m= 
o-= 6 C= ‘ Vi ee Ge == * oo te 
® @ ] bay ge ‘5 ~ > 
- ‘ . y , s 2 Pa om ve 

a 1 wm 48 (6G Mi Gs *evh 

‘ a | ) en oe 


oa 
- a : > on ; * iA a0 et ae vy 
5sea<2chk We Am Oe eA 
e > © be te be 6-9 OG ES CT : 
_ oe : 
ut ¢€ ot g@e or ie 4 17o <= & tu w 
° Oey 2) > le de —- 
~ ee? e& se 9'° vs. 4 
7 >) (tear oan | 
mnreooseee a = 


a 
_ 


- 
¢ OO ews WH Vv : We é 
“6 =e I - . ‘ . é ) 
a 
- & 
i 5+ Gee o ely « eta 
ad =» (07 -¢6 ee? Po, 4 
' aay = 
i. se =z 


0 
a 
& 


7 


-_ a * ag ,e@ ~ 
>  # f . % e » «& e an 

¢ ’ ? q 
a irs ‘a in ae) 


" 
, os Ny 9) O¢r~ & OG 
Pm yvie% Tr 
° ® = 
4 i (Als ny ae 
a ‘ @ 9) =o! Gs = 
tw? 7 3° i FA ts Fg Al 
> > - 
24400ene4o. A500 
ee ee a VE MA SB, 1 Sr Ve 
= aoe | =e > >= aes (es 
- tan ¢ 


UA ts S&S 2 OPP - “aree © 
| 


“aS - & pp aa 


7 
» 
® _ 


w 


fe 


Ost 
Low 


€6v 0 


OLem0 
SE¢cu0 


(a: 


esl 
oa 


SoU 


vv9 0 
Sem 


Cu 


Wey: 
£81 


pop 0 


Clicau 
2.0 tae) 


Ocl 


bot 
eye 


SESE 1) 


Ua om0 
Sys 78) 


801 


Oat 
ae! 


,3uU0d 


0v9 0 


GLO my 
LE I80 


q 34eq 


G2 oe iL Ga aN IN 


Li eee mee 
hele WAR (0 FEC 


sLaWee Loe 


™ 


Gat 


157 


8ST 
Gal 
mle 
Sve 
Lt¢ 
Else 
ook 
SIT 
Oa 
SOT 
86 

G8 

26 

Cel 
aust 
Vall 
cat 
6vl 
Sel 
AEM 
207 
Tot 
16 

vl 

Be 


JASW 


T= 
yz0 


vLcev 
SSE 0 
Oe 0 
O6T 0 
Bele 0 
p02 0 
6ve 0 
So hSea0 
ES 0 
S09 0 
shen, 
Ss () 
Boe} 
euro 
Grane 10 
LoS C 
6 ae 
oie 
score at 
faite! 18 
Sie 1 
£020 
89L°0 
660. 0 
650 0 

J 


SAYMN 
owsun 


oan com” 
CUA Oe 0) 
Vv oe ee 
Oral ee ie 2 
96 e391 0 
EE Vas 
WG eee 10 
69 Saag a0) 
6L 069°0 
BZ SOs 0 
vl pOL°0 
vl ae) 
76 Sie (8 
(uate. 
Ora ee re) 
97l vSv'0 
I feu 10 
HSL eee 18 
ewe alensn 16) 
Gros emer Orrin) 
jelly © WARNS 
Ud tae U CoeU 
Gat CEs 7a 
OS Pacey 0 
Chee hyo <0 
J4ASW J 
-TT=4AVMN 
Caen 
AdYy}OOWS 
QO) Lahti 


Oo cia) 
Co eG 9S 30 
Volpe cae) 
Ge MOE 16) 
OM = (a SISF 6 
6L ae eat) 
8S Spi 1) 
Ge, O7L0 
Lied Cee 0 
Sie ieyer 10) 
Tee Sach 0) 
ee (S/S a8) 
96 8Lv°0 
Gi Onecan.) 
Wo tenis Ae 
eek ashe 1G) 
Gai cca 
Orcs ieee cere 
oy Ova) 
ons en oy an aoe 
OT ear ce 
Ov vcs 0 
Oo eee) 
Gol 297 C 
coe 655.0 
JASW J 
9=4AVMN 
€ “ON 
AIYZOOWS 

= JAVMC 

= €WNN 


LN Neves 10 
Sole 409.0 
Gu laeec 0 oa) 
Ol eene Com 
Tee eve ok 0) 
Ss ale 
3 Dea 
sil GEO 
LCeaee sagan 
Ue eae) 
SC) CO mau 
69 evish 
EL es) 
Stl soo cO 
Bre | ASAT 18 
Cxes Ay 18 
Di ey Sa) 
ae ails 18 
WAN tens ier 38) 
Pet 7S S50 
See SUE @) 
Mes SESE 
Vi ncog. VU 
GoD ylece 
OG a0b9=0 
4JASW J 
9=3AVMN 
ee aN 
AIUPOOWS 


495 e1eg UYauaNOSs BuLspn szyNnsay 


ene Sse 
6 oo G0 
PME sy 
Reba 1elere) 8 
26 esvie 18 
09 COG 
Be C930 
Tp E50" 0 
ee c6L°0 
os chek 
UZ o Lan 
v3 p65 °0 
86 cee 
OSL Tod =O 
Sey Skars -710) 
OO lame Ounces 0. 
Mepe  G SCS 
Oo aac O cane 
pele NSPS, 19) 
CS ce. 0 
Orel OL, 0 
Se 876 0) 
pk = vey 0) 
Oey Vac 
Cee Sac 
JISW 4 
9=4AVMN 

tT °ON 
ADYPOOWS 


9 TeeG ce caw) 
De toe G0 0) 
Ora elias 0 
CVC. V0ce 0 
Oc Cem. aed 
GU Cae ie ced 
SUI ane a0 
OL 957-0 
eee ASW 
OUT 2 eo 
96 spec) 1 
18 039-0 
88 659° 0 
CO eect) 
Sb celeys 10 
CC ee U, Se0 
stake SS 0) 
O97 San Vie ge 
say h Se ie 
Chae aimee (0) 
c0T 08970 
€8 SCRE, 40) 
99 Tvs °0 
IOs 9 LOO 
62 856 0 
JASW J 
9=4AVMN 
YZOOWS Uf 
“TW 9L9gel 


(a) WQS 190) Cag (Ue) Wein Tee) (en) 


she 


4,3 
Be ies, 8 
Shi 


a é- - oe a 
om, Qe. pee z : . . ; vcr 2td,.8 
et et 3G eS 5G 2¢ ize, 
Gee 8f8 | = ap ti. Gea 
" Br? y 3 : - e° /GekeF 
. nga .@ < 4 c3 9 
_ oa es be Te 
=—s . J “a2 es, 
eo5 ¢ ; tak ita 
E¢< = : : nl “s ra). _ 
a . 


¢ ‘ 
oc! a : 7 : ea d£ a | er 


+e 


ize) 


vol 


661 
SCG 
© Ga 
Ost 


O87 0 


Oa Can 
col 0 
ef¢2 0 
Can 


MTA 


90¢ 
96T 
89T 
Gi |) 


eye) 


Hane 10) 
88S 0 
vv9 0 
ESTO 


Jer 


vO0¢ 
661 
Was 
O9t 


Sass) 1G) 


ree) 
eager il) 
Se eGeig) 
CAVED (0) 


oot 


ilk 
Care 
pst 
BOT 


ae 


bOJa0 
inresse 0 
0790 
679 0 


ae 


00¢ 
6L1 
oa 
€evl 


LoS 30 


Sree 1) 
Caen 
CUS a 
co9 U 


005 


2) IW) 


cvl 


c8l 
00¢ 
98T 
vst 


c87 0 


cue 0 
ors ta 28) 
T6¢2 0 
cise 0 


NWSW 


8é 
Le 
92 
Ge 


“TW 9PLgel 


SS oh 
es ues 


wt mt if 


2 is 
si ate uv aae isi 


aa 
af 


_ 
> 


i“ 


ik : 
7 


159 


Susy. ens Te Ad) 
Came Vn) 
SIO TSC 
CB) topes (G 
Doe Coe 0 
a (ial Tea, 
EG. 709.0 
Oi Pee 
Pie wae hee) 
Greece nO 
0 Gan 0 0 
Lem 030.0 
OR SP ew) 
Ope 0¢cL 0 
Se GATE 18) 
ahi Rae 
o Ve Oran) 
£76090 
Sage 8 
DUM LO Ce U 
(AS werisye 8 
8S Zvi 0 
ep tee 10) 
8 CYS (a 
655 759 "0 
4JISW A 
T=4AVMN 

€ ON 
AIYZOOWS 

eB ZAed 


oe SS a 
69 8600 
69 BOG 
oe) aes 1) 
69 67S 0 
v9 lee AW 
19 ee 6 
LS 7:6 S540 
0S 5090 
Ov GLO) 
OV 939.20 
OV See 
60 Eee 18 
i (pene) 40) 
She) 9755.0 
GL 16720 
OL Ge iG 
6% epee 13) 
Gas ies) 
tS CooaG 
09 CO aU) 
Ge, Lao 
LL (yaey 1 
BL Bes 0) 
vl Gee 
4JASW J 
T=4AWMN 
cae oh 
Aayuyoous 
8 = 
OV cs = 


66 62£°0 Oe eet S81 200-0 
Gomer lY <0 go ei G0 Esl Svo'O- 
Si 28710 ome G20 I8l 8t0°0- 
89 8595970 cameo. 681 £€80°0- 
9 €09°0 6S v69°0 foe CO m0 
Beh EE SGenoes.0 L191 v£0'0 
GC mic oO GE 10 Bo Ge 
6Guec19. 0 Ze NG C0: oe 7s0 
SG" 96570 GCurOl Ge0 Bie RE 
9p 959°0 oD ine 8 Ce Ga 
ii RES en cOn0 88 +0S‘0 
ie BREW €p 06S°0 Me SE OG 
aie GI Go eaicc 6 ie eee 
LG) 7950 v9 90 (ae a 
eG 9 SESE G9 +8r'0 Ce leo ca 0 
me  eegra 9 894°0 pet 102-0 
29 = ¥65°0 65 6°60 1 ee Can 
6G mero GeO Cia Ge S01 606-0 
Goues09 "0 Dip PW €6 20-0 
2 MG Sp Sao e ek Cow w 
Omen oe ey Pee I8  89S°0 
we He 26° G9/°0 69 299°0 
ie GAO” 655 a cL/20 Ov o0€8°0 
Zomees ea Coan Gi — EO 
Cou a0 6y ©. 9690 0 000°T 
FERS J 74Sw J 33SW J 
T=JAVMN 
I “ON L=3AVMN 
ADUPZOOWS YUPOOWS Uf) JgQUdTSisSvod 
JAW 
qWfiN aS eyeg usayzNOS BHuisp szLNsay 


98 G77 0 ve 
6L Simie OP Rete 
€8 GeV Occ 
SL E80 OF Sic 
08 L67'0 O¢ 
oe ee wy il 
Si ACE 1. esi 
52 sry 8p Ce 
ey see 10) Sel 
8P Seer) Si 
0S 62950 sit 
ee Gc UmOmec | 
Us ose We el 
TZ Lae 19, AE 
aq 6Co Oe Ot 


si 667°0 6 
89 Sis) Oy 58 
she Aas (ee aL 
OS sispe) eee) 
i epee OPS 
2g OE 0 aaa 
8S Seepeth 1S 
OZ LG FO Dare 
vl 66-950 aaee 
oe G9S°0 0 
JISW J 


ABotozpeur gy 


ecV Ware et 


Sera wes Be 7 
z=. ow > 
> 


7 on 
» | a 4 = , 
>a om em ¢ ~~» ar 
. a aa ee ci 7s 
. “ 
> 4 % 
® 
' 
a 
i, a: 
a 
J t— be 
j > 
) ® -— @s ie . 1“ 
' L - . | 7 —<_ > - 4% 
i , ° i] 5 ‘ et 7 tit 
’ =» * \, lt § ue 
S ; ; + eS © 
' 
et 
a | 
7 7 
4 - 
= 
> . 2 eS = 7 
ie La 
ty : ; 29 G2 9 Am FF 
om » pea 7 = a 
=> ee ‘oe «we j re 


_ a) (iwss hed 
oe 


i a a 
, 


160 


bs 


OL 
OL 
OL 
6L 


O92 U 


since 10) 
oloe0 
Se 
6S 0 


EY 


09 
19 
6% 
OL 


ee 0 


COGN 
Giencl 10 
eat) 
es 0 


he 


b6 
06 
L8 
96 


O83 Gae0 


9€7'0 
Slt’ O 
Shy 18) 
T8€ °O 


6S 


Oe 
iE 
GL 
08 


See 0) 


ewe 18 
00S °0 
Scie 
Hashine 1 


cel 


CelaL 
9cT 
Gal 
OL 


eT UO 


Aas 10 


Sau) 
Steed 
Shee 10 
ecUr0 


Teh The Wes 


89 


hye 
GZ 
vl 
08 


“OW 


19$°9 NV4W 
SIS;0 82 
cOVs0 ee 
OC a0 aoe 
£9v°0 Se 
cic ea 


161 


L6 Ltr oO 
Es Teas 0 
Ee Sias 18) 
89 065 0 
ja Loa 0 
v9 OES) 
89 Geo 0 
84 oe) 
ae Oa) 
St JO 50 
cv geo 0 
Tv et 0 
85 Soa) 
je, aan 1) 
88 857 0 
96 St 0 
66 807° 0 
86 OSv'O 
PUieecor 0 
76 Eig 
68 [eeer 15) 
98 esis 8 
v8 ON 0 
tl 719-0 
JASW J 
v=4AVMN 
T “ON 
49YUPOOUS 
g 


18 Goya 
lime eee a8 
Eg GLa 0 
Eo 865 0 
Lg Bo ot) 
OL Pee aon) 
ie eg) 
6% Siw) 
6 CHOSE 18 
Ev BU 1s) 
BE exons 1 
Tv COLO 
8Y pv9 0 
SP VS oe0 
95 ELGg 0 
she SHAS, 8) 
a UV S)) 
ae OSC 
6S Gin ee 1) 
oe fe lhe) (8 
v9 189 0 
og Rol AU 
09 Gee) 
ae) e790 

JISW J 
€=4AVMN 
€ “ON 
ADYPOOWS 

8 

t4ed Ove 


i} 


69 gio 0 
OG chips Ae 
Le Sieh Ne 
ao Comm 
£9 655-0 
Gel c8t 0 
hue Se) 
aye 02720 
89 Cay 
89 OS 7 e0 
che) cv 0 
cL 3o7) 0 
fife eo red 
SL 6S?) 0 
6L Sve 
Ge v0S 0 
69 elS 0 
8S Wwe) 
oS eee) Ee 
i ose 10) 
v9 Tega T8, 
Lo col U 
v9 €89 0 
oo Sie 1) 
JASW J 
€=4AVMN 
YZOOWS Uf) 
SE 
GW 


£8 vou 0 
SL ECO 
Gye v7s 0 
CL OLS sO 
ee IES ee 18 
Oz £9520 
oo Tee Gea 
Ga Sips Ae 
iy LOL 0 
Ce aie 18 
Ee Cae 
Ge voL 0 
ov 889 0 
97 COI nU 
oe as) 16) 
OS Oro ae) 
Lv 179 °0 
LY eps) 18) 
0S Peo 0 
Lo QeoF 0 
Lg Geo 
Sie Choe. 
ee boos 
89 EIS 0 
4iSW J 
c=4AVMN 

€ "ON 

Ady ZOOWS 


69 060 °0 
eS HESS AS 
Lo abs 1) 
OL COSa0 
L9 tS 0 
69 S3v 0 
79 E OeSe 18 
SHS Teveieia’ 
cv O29 =O 
LY iiek 0) 
Sv Oloyse 10 
ee) 67S 0 
GL prea 38 
SL onon ame) 
98 O8€ 0 
88 l0v 0 
EL 6570 
OL bse ae 
TZ 0vS 0 
NeW jh hes) NG, 
08 209.0 
68 Tao ae) 
98 See 0) 
Lae v7s 0 
4ASW J 
C=4IAVMN 
¢ ON 
ABUPOOWS 


VIAGme Ted UtdUINOSBOULS mise [nso 


vl Va OS Eo 
69 O85 0» 2é¢ 
Gee, Doo O ke 
oo (el 0h aes Oe Or 
09 609s 0amo. 
69 sag Oe rei! 
OL Gsp'o LI 
89 LV Oo. 
oS S09 0a 
a Vio oe. eee 
bs Voc Own 
og Cre ae 
G9 eeip th RAS 
8, Shenae thy 20) 
sh) O8U50 2.6 
6% LoS O ES 
97 Oo Oe 
ve BSR IO OS) 
ibe Ub Ome 
ore O¢3) 07 
one Leo 0) e 
Sv Gio 0) -¢ 
0S OS 207 
Gv Ey ORO 
4aSwW J 
c=4AVMN 
UPOOWS Uf} 
Ge We og ed 


SS 
ese tae ‘sts 
. _ _ - : 
an a : ™~ ; = 
a Ateseexd | gad toone 
> eagvaye . i“, oe 
a, Se “. tebgee? 


ei: 
ae" 


. 


“4 


a : ne 
%Fe ¥ *s 
— ee» J 7 
ee *o9 “Sp 
a oP.5 <5 
23 . “hE 

5 ae 
> 


Fie 688 (2) Op Ei “we o's 
acs 
a 


162 


wae ae) 


pee 0 
UO Lea 0 
New 
Roos 
Love 0 


OL 


83 
68 
18 
G8 
(byh 


Ova 0 


Soe 
cues)» 19 
sete 8 
Boa 
ieee 0 


OOT 


Bos 0 


vey 0 
86E°0 
Cy VaeU 
Woks 10) 
LOCO 


69 


is 
v9 
oy) 
Sue 
il 


Ses 0) 


soe @ 
69 a00 
Cam 
lov 0 
067 °0 


,3u09 


q y4ed 


ue) 


GL 
BL 
AVE 
v8 
G8 


G6S°0 NVAW 
609°0 82 
ee be Ie 
Gens, 0) eke 
AASV 0) SES 
Ogun Uavc 


"SW PLGeL 


Gs) 


€ “ON 
JayZOOUS 


L8 Sa 1 
EE [ere 
ck jase 8) 
pe cneken 0) 
ele Cae 8 
OL oS 18 
v9 Susie) 40) 
eo Chesee Ae) 
cv 2UL 0 
Ge SLE GD 
Ste eye 18) 
Ge SI 0 
90 98 980 
9” 7690 
Gs LE9° 9 
0S es 1) 
Lv 0v9 0 
LY pe 8) 
a4 eek 
oo aes @ 
us PLO e0 
Nil oan) 
a [Eee 
ne yee 8) 
JASW J 
TL=SAVMC 

¢ =IAVMN 

€ °ON 
AdUPOOWS 


a adie ine = 


8L c8t 0 
69 Ly Sa 
ee) O75 0 
89 eas 1), 
Oz Gat 18 
69 GeO) 
oo OS. 0 
0S TS 9e0 
Ov CEO 
ce eet 
ve vvl 0 
on shea 0) 
ts Ly oa 
om ue 0 
8S Gee ae 
GS eles) 16) 
18S £6570 
OS eae: (0) 
oo Siders) 10 
89 Syeer 40) 
tte Gu, We 
03 Ly oat 
Sif ches, 10) 
ng CU Sag O 
JASW J 
GS=4AVMC 
C=4IAVMN 
one oh 
ADYUPOOWS 
aWAN 


d Oe Seo 
89 isewie| 10) 
ee 009" 0 
Je) oe 0 
ve OU Sea 
vl Lov 0 
GL SL 0 
Ga fea 0 
aie HES 16) 
St Sai) <9) 
ay Teoma 
6€ ale AU) 
iS 90.950 
as 8E9°O 
OL Wecack 8) 
Git oe 8) 
LL cca 0 
69 case 
OL SLS 0 
IZ reo ma 
ee BO 
ie opie 718 
8S OPA 40 
SPs 669°0 
J4SW J 

S=4AUMC 
G=3AVMN 

ce ON 
AGUZOOWS 


18 6907 °0 
LL Lovo 
bye O50 
Sy Ee 8 
iki Gi a0 
8L Gey) 8 
Ore vov 0 
to ios 10 
60 OL oe0 
cv G39 0 
Ov Sogeo 
6€ Cee 0) 
0S ess) AU 
vS EVO 
oes) USS °0 
ere) eo 
en) ba 10 
09 eens. 1 
o9 Sees 18) 
Ke eek 1G) 
6S LOU 
6% LV 
vs 602L°0 
TS S69. 0 
4J3SW J 
S=4AVMC 
v=4AVMN 

€ “ON 
ADYZOOWS 


29S e}eq Usau NOS BHuLispn sa[Nsay 


6L 9bb'O0 Ez 
08 Mmeavecmecd 
nS ae 2 
28 «28°02 
6£ v0S'0 61 
9f €0G°0 8 
pL €05°0 LI 
19 €99°0 91 
€g so9°0 SI 
pS 8090 oI 
[Gere (og meet 
29 yd hoa 
LL 8dv'0 IT 
08 9°00 OT 
06 6€°0 6 
gg €€p'0 8 
9£ g9p'0 Z 
89 Orso 9 
99 6/5°0 &§ 
IZ g09°0 ¥ 
pL 2970 € 
€g8 92970 2 
28 «6p9S'0 1 
ae ie 
SW 4 
8=3AVNC 
p=4AVMN 
2 ‘ON 
ADYZOOWS 
“2°W ALdeL 


ose ~~ ws ~~ > 9 , 
i6 e 4 e/ O° Here , 
J 
> (@ ( ry on Ce 4 ‘4 
— is ne <> ¢ 0 ve 
' 2 x 2 UP ns ue - 
, >< — >< > & <p iF 
7 beli | g Y BD ee Dim @ vt Cc, & wo 
: ‘ ; : tho Sew POP HAF Che Oy SO tp) = 8 
»@ ry oa vw —_ & 
= » re Viae @ o Co oe 
e | ‘s = Vy Oy we nl Ole 
- ) 1 ow B ? aG 
i] 
* 
Se , j “4 (14.0 Bee 4 
f 4 t 4a ‘9G “ber ep OL e-* we? ER) Cie CO) VED 
4 tae Cw 4 @ *@ 
( ‘ >) i. > rr we ow oy 8 obs @ 
, ' ' Hj q) >) 417% ¢) Ge @ if 7 
* ’ ry @ 
> : 7 ¢ 4 => > 6 aw 2. art 5 . 
os = ap =_—A,) A P' rome > ° 
Ge “4 - =) $ a) ap Ge GO Ty 
As 9 ile @ j - (Ged > ort 1D Pm 
icc’ @ © n i (hm 1) os @ —= = «& Ls 
| ad : ad w SY . yw sw 4.7 
2 9 
P47 GO J Ga < » oy 2 GO  & € es &y 
©) #11 


SELES Shek Wiliam Aw NES Bees 


i, 
baie 


- » 
po * 
(akc im ee 
Qe ey 
zeae 
= ae 

i) 

' 

a 

wa we 

6 
Eero 

ee” 
+s B 
om 

- 
am 
so ve 


ogi 


#OKh = Sen 
y 


Hesse sapresiu-perse 


— 
oe 


= 
> 


164 


aie) 


76 
16 
68 
66 


O01 


(ee) 


Hiren 518) 
aU ve) 0 
6077 0 
698-0 
68€°0 


oo 


c8 
(ey 
Sl 
L8 
L8 


TES, 18 


v67 0 
vet 0 
Breer 
bare 39) 
Se 10) 


69 


SOT 
EOT 
16 
€6 
L8 


See 


627 0 
vivo 
8870 
OLV 0 
VU eet) 


L9 £9S°0 OL 
OOo e2Or = 0 16 
G6 ae 10) 06 
fants} 1 on) 88 
NS) Sasa 18 06 
06 097° 0 L8 
2095 op aera 


o0a80 


Str 0 
600° 0 
vor 'o 
607 0 
6€v 0 


NV AW 


8c 
ve 
9¢ 
G¢ 
ve 


"?°W 9LGel 


sake 


= 
anes 


, ary 
Pe 
w mB mits 


APPENDIX N 


Actual Forecast Results 
(Northern Data Set) 
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Results Using Northern Data Set NUMB = 240. 
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